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Close  to  the  onset  of  self>excited  fluid  oscillations  the  generic  complex 
Glnzburg'Landau  is  proposed  as  the  lowest  order  model  for  the  "plant”.  Its  linear 
part  which  provides  the  stability  boundaries  is  derived  from  first  principles  for 
both  doubly- infinite  and  semi-inlinlte  flow  domains.  Concentrating  on  a  single 
global  mode,  the  model  is  further  simplified  to  the  Stuart-Landau  equation.  For 
this  latter  model  a  methodology  is  developed  for  the  design  or  single-input 
single-output  controllers.  The  so  designed  controllers  have  been  implemented  on  a 
self-excited,  heated  two-dimensional  Jet  with  one  hot  wire  as  sensor  and  an 
acoustic  .peaker  as  actuator,  and  are  shown  to  be  effective  within  their 
limitations  in  suppressing  or  enhancing  limit-cycle  oscillations.  Finally,  the 
effect  of  of  a  controller  designed  to  suppress  tlie  most  unstable  global  mode  on 
other  modes  is  investigated  experimentally  in  the  wake  of  a  cylinder  at  low 
Reynolds  number,  where  an  encouraging  semi-quantitative  correspondence  to  the 
Ginzburg -Landau  model  is  found. 
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1.  THE  DERIVATION  OF  THE  LINEARIZED  GIKZBDRG-LANDAU  (GL)  EQUATION 
FROM  FIRST  PRINCIPI^S 

1.1.  Related  Publications  emd  PresentatioT:^ 

MONKSWITZ,  PETER  A.,  HUERRE,  PATRICK  &  CHOMAZ,  JEAN-MARC,  Global  linear 
stability  analysis  of  weakly  nonparallel  shear  flows.  Submitted  to  J.  Fluid 
Mech. 


1.2.  Introduction  and  Basic  Equations 

This  part  of  the  report  addresses  the  extension  of  the  concept  of 
absolute  instability  to  nonparallel  flows.  As  in  the  case  of  parallel  flows, 
in  which  an  initial  impulsive  excitation  leads  to  the  ultimate  dominance  of 
the  most  amplified  Fourier  or  normal  mode  with  zero  group  velocity  (Briggs, 
1964,  Bers,  1983,  euc.),  we  attempt  to  answer  the  question  of  the  asymptotic 
(in  time)  Impulse  response  of  a  nonparallel  flow  and  we  will  call  the 
equivalent  of  the  "absolute"  normal  mode  a  linear  "global"  mode.  The  latter  is 
simply  a  time -harmonic  solution  of  the  homogeneous  linearized  disturbance 
equations  with  homogeneous  boundary  conditions  in  space.  Such  solutions  can  in 
general  be  obtained  only  numerically,  especially  if  the  basic  flow  is  strongly 
nonparallel.  Examples  of  such  computations  in  shear  flows  have  been  published 
by  Zebib  (1987),  Hannemann  &  Cartel  (1989)  and  many  others. 

If  the  the  mean  flow  is  weakly  nonparallel,  i.e.  evolves  slowly  on  the 
scale  of  a  typical  instability  wave  length,  global  modes  become  accessible  to 
WKBJ-type  analyses.  This  involves  the  step  from  the  "slowly  diverging" 
approach  of  Bouthier  (1972)  and  Crighton  &  Caster  (1976),  who  treat  the 
spatial  evolution  of  a  forced  wave  in  an  inhomogeneous  mediiun  ^the  signaling 
problem) ,  to  the  problem  of  finding  the  unforced  global  modes  vmare  the 
streamwise  direction  also  becomes  an  "eigenvalue  direction" . 

In  the  following  we  treat  the  case  of  infinte  or  semi -infinite  shear 
flows  that  contain  regions  of  both  absolute  and  convective  instability  by  r 
WKBJ  analysis.  In  particular,  the  cross-stream  structure  the  global  noJes 
is  taken  into  account  and  an  explicit  discussion  of  the  connection  hetveei  the 
properties  of  global  modes  and  local  absolute  and  convective  instability  l'> 
given.  This  connection,  which  is  supported  by  several  examples  (Moijkswitz 
1990),  is  put  on  firmer  ground.  We  note,  however,  tha»  the  analysis  admltf/ 
only  "local  feedback"  by  vortioity  waves  to  drive  global  modes.  It  is 
therefore  restricted  to  cases  where  long-range  feedback  'a.  negligible  auc‘. 
important  problems  such  as  edge  tones  are  not  addresfati,  £or  which  the  primary 
driver  of  the  instability  is  the  acoustic  feedback  fro*c  a  downsti-am  point  of 
intense  fluid-surface  interaction  to  a  trailing  edge. 

To  study  two-dimensional  instability  waves  in  c  spatially  inhom  ganeous , 
incompressible  medium,  we  start  from  the  equation  for  the  z-vorticity  -V^  , 
where  ’J'  is  the  total  stream  function. 

[3^  +  (ay»)a,  -  (a^'i')3y]  (2.1) 

Next,  'P  is  decomposed  into  time -independent  mean  flow  ip  and  a  small 
disturbance  i}>' 
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*  -  ^(X,y)  +  rf>'  ;  X  -  cx  .  (2.2) 

At  this  point  we  assume  that,  in  the  terminology  of  the  method  of  multiple 
scales  (see  e.g.  Bender  &  Orszag,  1978),  the  mean  flow  ^  depends  only  on  the 
"slow"  coordinate  X-ex.  The  parameter 

«  -  A,y^{r"(x)(d5/dx])^^  «  1  (2.3) 

characterizes  the  degree  of  spatial  inhomogeneity  of  the  basic  flow  by 
providing  a  measure  of  the  small  change  of  the  typical  cross-stream  length 
scale  5(x)  over  one  typical  instability  wavelength  A  .  In  the  absence  of 
body  forces,  the  assumption  of  a  slow  evolution  of  the  mean  flow  immediately 
restricts  viscous  effects  to  0(e),  or  in  terms  of  the  Reynolds  number  X 

X  -  e'^R  ;  R-O(l)  .  (2.4) 

Introduction  of  the  decomposition  (2.2)  and  of  (2.4)  into  the  governing 
equation  (2.1)  first  leads  to  the  boundary  layer  equation  for  the  basic  flow 

(3y^)Ojj3j^)  -  (3jj^)(3j^)  -  R-'  (aj^)  (2.5a) 

b(X,y)  -  ay\i  ;  V(X,y)  -  -  -e'^  (2.5b) 

Linearizing  around  the  basic  flow  and  keeping  only  terms  up  to  0(c)  then 
yields  the  following  equation  for  the  small  disturbance 

[(a,  +  u  a,)?"  -  (a5u)aj  r  + 

e[v  8  +  (3^ayU)aj,  -  R-^V^V^]  r  +  0(€^\yl>'\,\r\^)  -  S(x,y,t)  , 

where  a^^  and  the  Laplaclan  have  not  yet  been  split  into  fast  and  slow 
parts.  The  source  S  has  been  added  for  the  study  of  the  impulse  response  in 
the  next  section,  but  the  ultimate  aim  of  the  paper  is  the  search  for  global 
modes,  i.e.  homogeneous  solutions  of  (2.6)  with  homogeneous  boundary 
conditions  in  space. 


1.3.  The  WKBJ  Approximation  for  the  Green  Function  emd  its 
Breeikdown 

Following  Bouthier  (1972),  Crighton  and  Caster  (1976)  and  others,  we  use 
the  WKBJ  approximation  up  to  the  level  of  "physical  optics"  (see  e.g.  Bender  & 
Orszag,  1978)  to  describe  the  evolution  of  a  pulse  on  the  weakly  nonparallel 
basic  flow  and  identify  the  locations  of  its  breakdown  which  are  "turning 
points"  of  the  problem  (see  section  10  of  Bender  &  Orszag,  1978). 

To  obtain  the  Green  function  G  of  (2.6),  the  source  at  x-x*  and  ^-y*  is 
specified  as 

S(x,y,t)  -  {5(x-x‘)  +  i[}r(x-x*)]’M5(y-y*)5(t)  .  (3.1) 

This  form  explicitly  accounts  for  the  non-analyticity  of  G  on  the  imaginary 
axis  of  the  wavenumber  plane  when  the  lateral  extent  of  the  flow  domain  is 
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infinite,  as  discussed  in  detail  by  Huerre  &  Monkewitz  (1985).  Next,  we  take 
the  Fourier  transform  of  (2.6)  in  time  according  to 

G(x,y,t)  -  (2jr)"^  J  G(x,y,w)  exp(-i£.>t)  dw  ,  (3.2) 

L 

where  the  contour  L  is  taken  parallel  to  the  real  a>-axis  and  above  all 
singularities  in  order  to  obtain  a  causal  solution  (Briggs,  1964).  Finally, 
all  x-derivatives  in  the  disturbance  equation  (2.6)  are  transformed  according 
to  the  chain  rule 

^  •  <3'3) 

keeping  in  mind  that  3^  and  do  not  commute.  Concentrating  on  a  flow  domain 
of  doubly- infinite  screamwise  extent,  G  is  required  to  vanish  at  up-  and 
downstream  infinity  as  well  as  on  the  lateral  boundaries  jyl*^  .  Hence, 
following  section  10.3  of  Bender  &  Orszag  (1978),  the  WKBJ  approximation  for 
G  can  be  written  in  the  standard  form  away  from  the  source  at  x-x* 

A  A  A  X 

-  {G*(X,y)  +  <Gf(X,y)  +  0(e^)}  exp(ir^J‘  k*(X>  ;u>)  dX' )  ,  (3.4) 

X* 


where  e  plays  the  role  of  the  WKBJ -parameter.  The  superscripts  "+"  and 
denote  the  approximation  downstream  and  upstream  of  the  source,  respectively. 
The  k*(X'  ;ii))  are  the  corresponding  local  wavenumbers  in  the  upper  and  lower 
half  k-plane,  respectively,  as  shown  on  figure  2b  of  Huerre  &  Monkewitz 
(1985).  For  simplicity  we  assume  here  that  there  is  only  a  single  pair  of 
eigenvalues  k*. 

Introducing  the  WKBJ-Ansatz  (3.4)  into  the  Fourier- transformed  equation 
(2.6)  the  stability  problem  reduces,  at  leading  order  in  e ,  to  a  streamwise 
succession  of  locally  parallel  problems  which  are  governed  by  the  homogeneous 
Rayleigh  equation.  Its  solution  depends  only  parametrically  on  X  through  the 
shape  of  the  local  mean  velocity  pjjofile  U(X,y)  and  yields  k*(X;a>)  as  well  as 
the  local  transverse  structure  of  G^  up  to  an  unknown  amplitude  Ap(X).  The 
latter  describes  in  essence  the  "transmission"  of  the  insuability  wave  from 
one  locally  parallel  region  to  the  next  and  will  be  determined  at  the  next 
order. 

0(e°):  X(^*:k*,«,X)  -  0  :  |^J(|y|-«o;X)|  -  O  ,  (3.5a) 

with  G*(X,y)  -  A*(X)  ^^(y;X)  .  (3.5b) 

The  Rayleigh  operator  £,  which  includes  a  list  of  all  relevant  parameters  in 
the  argument,  is  defined  as  (see  also  appendix  A) 

i:(.:k,a>,X)  «  (kU(y:X)  -  mJOj  -  k^].  -  k[3jU(y;X)]‘  .  (3.6) 

At  linear  order  in  e ,  the  following  inhomogeneous  Rayleigh  equation  for 
G*  is  obtained,  where  the  derivatives  of  the  operator  I  are  defined  in 
appendix  A 
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0(£^): 


+  i  3jjk  /j^(^*;k*,«.X)  +  /^(^*;k*,a>.X)  }  . 


(3.7) 


It  will  prove  useful  to  replace  d^4>\  the  expression  (C.4)  of  appendix  C 
where  the  functions  etc.  are  defined  by  (C.5).  In  order  to  avoid  secular 
teinns  in  the  solution  of  (3.7),  the  solvability  conditii.a  (3.8)  has  to  be 
satisfied.  After  dividing  by  L^(^g;k*,«,X)  and  tising  (B.4)  for  the 

condition  reads 


a^u>*  a/*  -  -A*  (  iato*  +  (Va)^  a^^k*  +  u>^  a,^*)  .  0.8) 

a«W)  -  i  [L  (^*;k*,«.X)  -  l^(^^;k*.«.X)]  [L„(^*;k‘,a,.X))-^  ,  (3.9a) 

u>^(u>.x)  -  (2I^(^^;k*.w,X)  -  I^(^*;k*.«,X))  (L„(^5;k*.«.X) )-^  ,  (3.9b) 

«*^(«.x)  -  I^(^*„;k‘,a).X)  [L„(^*;k^a).X)]-^  .  (3.9c) 

The  L's  above  are  defined  by  (B.3)  in  appendix  B  and  represent  inner  products 
of  the  corresponding  £'s  with  the  solution  of  the  adjoint  Rayleigh  equation. 
Equation  (3.8)  is  now  readily  Integrated  to  yield 


A*(X)  -  a!*  exp 


r  A 

ij 

y* 


X  ( +  (Vp  ajjk‘  +  i  a^*]  (X'  ;w) 


ak«*(X'  ;w) 


dX' 


(3.10) 


This  brings  the  approximate  description  of  the  Green  function  away  from  the 
source  to  the  X^vel  o^  "physical  optics".  What  remains  to  be  done  is  the 
connection  of  and  gJ  across  the  source.  For  uhis  the  analysis  of 
Bender  &  Oj^szag  (.}.978)  has  to  be  generalized;  instead  of  patching  the  two 
solutions  G”  and  G*,  which  is  only  possible  for  ODE's,  we  use  in 
essence  the  approach  of  gurridgp  and  Weinberg  (1977)  in  which  the  connection 
is  achieved  by  matching  G~  and  gJ  to  the  local  parallel  Green  function 
at  X®.  The  latter  is  obtained  from  the  double  Fourier  transform  6  in  time 
t  and  space  x  at  X-X®  which  is  given  as  equation  (8)  in  Huerre  &  Honkewitz 
(1985).  Integration  cf  G  with  respect  to  k  along  the  contours  shown  on 
figure  2b  of  Huerre  &  Honkewitz  (1985)  yields,  upon  evaluating  the  residues  at 
k"*^  and  k~  respectively  for  any  frequency  on  the  contour  L  of  (3.2), 


2H(x-x®)  ^;(y;X®)  ^;(y®;X®)  exp[ik"(w:X®)x] 

G(x,y,u>;X®)  -  - - + 

[k^w:x”)U(x®,y®)  -  U)]  a^D[«.k^u:x“)] 

2H(x®-x)  ^Q(y:X®)  ^'(y'iX®)  exp[ik"(«:X®)x) 

+ 

[k-(w:x®)U(x*.y®)  -  «]  aj^D[u..k"(«:x®)] 


(3.11) 


In  this  expression  H  is  the  Heaviside  function  and  D  the  dispersion  relation 
associated  with  the  local  Rayleigh  equation  at  X®.  We  note  that  in  (3.11)  the 
contribution  from  the  integration  along  the  imaginary  k-axis,  which  is  a 
branch  cut  of  the  dispersion  relation,  has  been  omitted  under  the  assumption 
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that  the  long-time  behavior  is  dominated  by  ♦’he  discrete  spectrum.  Comparison 
of  (3.4),  (3.5)  and  (3.10)  with  (3.11)  finally  yields 

2  ^*(y':x') 

a;*  -  qp -  .  (3.12) 

[k*(a>;X*)U(X*,y')  -  «]  a^D[o),k^(a>:X') ] 

We  now  turn  to  the  discussion  of  the  long-time  behavior  of  G  which  runs 
analogous  to  the  discussion  of  absolute  and  convective  instability  in  the 
parallel  case.  The  basic  idea  is  that  the  leading- order  time;; asymptotic 
behavior  of  G  is  determined  by  the  uppermost  singularity  of  G  in  the 
(i>-plane,  i.e.  the  singularity  with  the  largest  temporal  growth  rate,  which 
becomes  "pinned"  as  the  tj- contour  is  lowered  .  In  the  parallel  case,  the 
singularities  in  the  co-plane,  which  correspond  to  zeroes  (o(k)  of  the 
dispersion  relation  D,  can  be  moved  by  deforming  the  Fourier- inversion  contour 
in  the  k-plane  until  the  latter  is  pinched  between  two  branches  k'''(co)  and 
k"(co)  (see  Briggs,  1964,  Bers,  1983,  etc.).  When  this  happens,  the  singularity 
to  ■<o(k  )  becomes  "pinned"  at  the  absolute  frequency  to^  which  corresponds  to 
t^e  saddle  point  k^  where  the  complex  group  velocity  is  zero  and  the  k- 
cortour  is  pinched. 

In  thg  weakly  nonparallel  case  we  can  argue  in  a  completely  analogous 
manner:  .'f  G(X,to)  becomes  singular  at  a  location  X‘  from  which  the  X- 
integratios'.  contour  cannot  be  moved,  the  corresponding  pole  to(X*')  becomes 
"pinned".  Again,  the  "pinned"  pole  with  the  largest  temporal  growth  rate 
determines  the  time -asymptotic  behavior  of  G.  From  (3.10)  it  is  seen  that  G, 
in  particular  Ap(X) ,  becomes  singular  at  the  zeroes  of  Such  points 

are  in  fact  turninf;  points  x‘  where  the  WKBJ- approximation  breaks  down,  which 
is  easily  seen  by  noting  the  correspondence  between  Ap(!0  and  [Q(X))*^^*  in 
standard  textbook  notation  (see  chapter  10  of  Bender  &  Orszag,  1978).  In  the 
following,  we  have  to  distinguish  between  two  types  of  pole  "pinning": 

The  first  is  the  exact  analogue  of  the  parallel  case  where  the  X- 
integration  contour  is  pinched  between  two  branches  X*(«|3j^fa^0)  on  which  the 
group  velocity  is  zero.  The  pinching  occurs  at  the  saddle  point,  or  second 
order  turning  point  X*"  where 

3jjW^(X‘)-0  ,  (3.13) 

as  shown  on  figure  lb.  This  is  the  case  discussed  in  detail  by  Chomaz,  Huerre 
&  Redekopp  (1991),  which  arises  when  the  absolute  growth  rate  w^j(X)  has  a 
maximum  within  the  flow  domain,  and  represents  a  generalization  of 
Pierrehumbert' s  (1984)  frequency  selection  criterion  to  the  complex  X-plane. 
The  assumption  that  a  maximum  of  w^iC^)  exists  within  the  flow  domain  Is  in 
fact  rather  weak  and  corresponds  to  the  existence  of  a  maximum  of  the  temporal 
amplification  rate  and  an  associated  saddle  point  of  w(k)  in  parallel  flows 
(see  Caster, 1968) .  In  the  following  we  make  the  same  assumption  and  consider, 
for  simplicity,  only  one  saddle  (3.13)  between  regions  of  stable  flow  far 
upstream  and  downstream. 

The  second  possibility  of  pole  "pinning"  arises  when  the  complex  group 
velocity  first  becomes  zero  at,  say,  the  flow  boundary  X“0.  In  this  case  the 
pole  w  (X‘')  is  "pinned"  because  the  X- integration  contour  has  to  start  or  end 
at  the  boundary,  as  shown  on  figure  Ic.  Hence  we  have 
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-  0  (3.14) 

and  in  this  case  corresponds  to  a  first  order  turning  point.  This  is  the 
situation  we  will  consider  for  the  semi> infinite  domain  where  the  upstream 
boundary  at  X-0  dominates  the  evolution  of  the  disturbance  while  we  assume 
that  the  flow  is  stable  far  downstream.  The  same  situation  has  already  been 
considered  in  a  model  by  Chomaz,  Huerre  &  Redekopp  (1988). 

In  both  cases  the  global  mode  frequency  is  given,  to  leading  order,  by 
w  (X^)  and  global  instability  is  determined  by  the  sign  of  the  leading  order 
global  growth  rate  u^^(X^).  Hence  it  appears  that  the  mean  flow  must  contain  a 
region  of  absolute  Instability  for  a  global  mode  to  become  temporally 
amplified.  This  is  confirmed  by  the  detailed  analysis,  which  also  yields  the 
next  approximation  of  the  global  node  frequency  beyond  w^(X^) .  The  latter  is 
obtained  as  eigenvalues  of  the  streanwise  two-point  boundary  value  problem  and 
are  determined  either  by  the  connection  of  the  upstream  and  downstream 
solutions  throu^  the  second  order  turning  point  (3.13),  or  by  the  connection 
of  the  downstream  solution  through  the  first  order  turning  point  (3.14)  to  the 
boundary . 


Xr 


Xr 


Figure  1.  (a)  Sketch  of  the  Fourier  inversion  contour  L  in  the  complex  w- 
plane,  (b)  The  WKBJ  integration  path  M  in  the  complex  X-plane  for  the  doubly- 
infinite  flow  domain  with  images  X'*^  and  X'  of  L  on  which  the  group  velocity  is 
zero,  (c)  Corresponding  path  N  and  X*  for  the  semi- infinite  case.  The  top  row 
(subscript.'?  1)  shows  the  situation  where  L  is  above  all  singularities  (O) , 
while  the  bottom  row  (subscripts  2)  shows  the  points  (0)  of  the  X-contour 
where  the  breakdown  of  the  WKBJ -approximation  can  no  longer  be  avoided  as  L  is 
lowered. 


I 


1.4.  The  Turning  Point  Region  for  the  Doubly-Infinite  Domain 
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In  this  section  the  ideas  of  Soward  &  Jones  (1983) ,  Huerre  &  Monkewitz 
(1990)  &  Chomaz,  Huerre  and  Redekopp  (1991)  on  global  inodes  in  a  doubly 
infinite  flow  domain  are  applied  to  shear  flows.  In  terms  of  application  the 
following  analysis  may  not  only  be  useful  In  truly  infinite  domains,  typically 
found  in  geophysical  shear  flows,  but  also  in  finite  flow  domains,  as  long  as 
the  boundaries  do  not  significantly  influence  the  flow  instability.  The  two- 
dimensional  bluff -body  wake  appears  to  be  a  case  in  point  as  suggested  by 
Monkewitz  (1988)  and  the  numerical  experiment  of  Triantafyllou  &  Kamiadakis 
(1990)  who  obtained  essentially  the  same  KArmAn  vortex  street  after 
replacing  the  cylinder  by  an  inflow  boundary  condition  downstream  of  the 
cylinder. 


At  the  turning  point  defined  by  (3.13)  the  first-order  equation  (3.8) 
for  A*  becomes  singular  and  one  has  to  bring  in  the  second  derivative 
9^.  Since  X*"  is  also  a  saddle  point  of  the  absolute  growth  rate  u^(X) , 
(Wq-w^)  behaves  like  (X-X*')^.  Hence  X^  is  a  second  order  turning  point°where 
the  second  derivative  of  A*  must  be  of  the  same  order  as  (X-X*’)^A*  (see 
Bender  &  Orszag,  1978).  This  leads  immediately  to  the  rescaling 

X  -  €‘^'^(X-X‘)  .  (4.1) 


In  the  inner  turning  point  region,  characterized  by  |X|:S0(1)  ,  the 
disturbance  streamfunction  is  expanded  accordingly: 

-  (*0  +  +  0(<®'2))(X,y)  X 

exp(i<'V(X-X‘)  -  iwgt)  , 

where  k^"k  (X*')  .  By  the  same  token  the  global  frequency  is  expanded 
around  w  (X')-<»»‘ 

,  (4.3) 


and  the  mean  flow  components  are  expanded  in  Taylor  series  around  X^.  In  terms 
of  the  inner  variable  (4.1)  they  are  given  by 


U(X,>)  -  U(X\y)  +  e^'^X(3jjU(X\y))  +  e  (x2/2)  (35u(X\y) )  +  , 

V(X,y)  -  V(X\y)  +  0(€^'2)  . 


Introducing  the  expansions  (4.2)  and  (4.4)  as  well  as  the  "slow" 
variable  (4.1)  (using  the  chain  rule  into  the  governing 

equation  (2.6)  yields,  at  leading  order  in  £,  the  local  Rayleigh  equation  at 
X*' 

0(£°):  /‘(^J)  -  0  :  |^‘(|y|-®:X)|  -  0  ,  (4.5) 

where  ^o(y)  is  the  local  eigenfunction  and  the  short-hand  notation  (see 
A.l),  together  with  the  corresponding  integral  (see  B.3),  is  defined  by 

L*^(‘) 


B  /(•;k‘,«^.X")  ,  (4.6a) 

B  L(.:k^,w^,X‘) 


(4.6b) 
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Hence  the  leading  order  solution  In  (A. 2)  Is  given  by 

*„(X,y)  -  Ao(X)  ^^(y)  .  (4.7) 

with  the  "free"  amplitude  A.(X)  to  be  determined  at  higher  order.  At  the 
next  order  0(e^^*)  a  solvability  condition  yields  only 

Wj  -  0  .  (4.8) 

Hence  we  have  to  proceed  to  0(e)  to  determine  A^  and  obtain  the  solvability 
condition 

3|Ao  (0^/2)  +  ix  ajA,  +  Ajz;^  -  -  x2(<4/2)]  -  0  . 

(4.9) 

with  the  abbreviations 


<4  -  [2LJ(4)  -  l4(^5)][Li(^J)]-^  .  (4.10a) 

4  -  (1^(4)  + 1^(4)  ■  4<4)(L^4r'  .  (A.iob) 

5«’^-1(lJ(4 . 1‘(4))[L‘(4]''  ,  (4.10c) 

4  -  (21^(4)  .  4(4)ilx>]''  .  (^-lod) 


where  the  ly,  etc.  are  defined  analogous  to  (4.6b).  The  resulting  amplitude 
equation  (4.9)  is  a  linearized  Ginzburg -Landau  equation  with  variable 
coefficients.  The  nomenclature  (4.10)  becomes  Immediately  transparent  when 
equation  (4.9)  Is  tranformed  to  the  spectral  domain  by 

ajj  -  i«*^'2(k-k‘)  .  (4.11) 

Using  (4.3)  and  (4.8)  this  yields  the  Taylor  series  representation  of  the 
dl.sperslon  relation  in  the  neighborhood  of  the  turning  point  X^ 

e  Wj  •*  *•  e  +  (4/2)  (k-k^)^ 

4(k-k^)(X.X*^)  +  (4/2)(X-X’^)2  ;^2) 

Hence  we  have  shown  that  the  dispersion  relation  (4.12)  which  had  been 
postulated  by  Huerre  &  Monkewltz  (1990)  and  Chomaz  et  al.  (1991)  is  generic  to 
the  turning  point  region  of  the  WKBJ  approximation  and  can  be  derived  in  a 
rational  fashion  from  the  governing  equations  under  rather  weak  assumptions. 
The  solution  of  (4.9)  is  now  easily  found  by  transforming  it  into  the  standard 
Hermite  equation.  Setting 

Ao(X)  -  exp((i/2)k^  5?)  0(0  with  $  -  (4a.^/w^)''*  X  ,  (4.13) 

(4.9)  becomes 

a^Q  +  a{[w2-5«V(i/2)4k^][4w^)-'^'=  -  CV4)  -  0  . 

(4.14) 

|q|  (?-+«)  “  0  . 


I 
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The  boundary  conditions  are  chosen  such  as  to  ensure  the  matching  to  the 
subdominant  WKBJ  solutions  both  upstream  and  downstream  of  X*"  and  restrict  the 
frequency  correction  to  a  set  of  discrete  eigenvalues  given  by  (4.15a) 
with  the  corresponding  global  eigenfunctions  given  by  (4.15b) 

-  exp[-fV^}  He^(0  .  (4.15b) 

where  tdie  He^(()  are  Hermite  polynomials  as  defined  by  Abramowitz  and  Stegun 
(1965).  We  reiterate  here  that,  according  to  the  discussion  leading  to  the 
frequency  selection  criterion  (3.13),  the  global  mode  (4.15)  represents  the 
asymptotic  solution  for  long  times.  The  recent  feat  of  Hunt  &  Crighton  (1991) 
who  determined  the  exact  Green  function  of  (4.9)  puts  us  into  the  unique 
position  of  verifying  this  statement  explicitly.  If  the  limit  t-^o  is  taken  in 
their  expression  (36)and  (44)  for  the  Green  function,  one  indeed  recovers  the 
most  unstable  global  mode  (4.15)  with  n-0.  The  higher  modes  are  more  stable 
since  the  imaginary  parts  of  and  are  negative;  corresponding  to  a 
high  wave  number  "cutoff"  and  to  stability  at  |X|-^  respectively.  At  this 
point  the  solution  of  (2.6)  in  the  inner  or  turning  point  region  is  complete. 
Its  matching  to  the  WKBJ- "tails"  presents  no  further  problems  and  is  described 
in  the  preprint  listed  in  section  1.1.  To  avoid  misunderstandings  it  is  worth 
pointing  out  here  that  the  notion  of  WKBJ- "tails"  does  not  in  any  way  imply 
that  the  amplitude  of  the  global  mode  should  peak  near  the  turning  point  which 
acts  as  the  "wave-maker"  for  the  entire  flow.  Depending  on  the  imaginary  part 
of  k^  and  the  downstream  evolution  of  Im(k'*')  the  wave  "leaking"  from  the 
"wave-maker"  region  can  experience  substantial  spatial  amplification. 


1.5.  The  Turning  Point  Region  for  the  Semi^Infinite  Flow  Domain 

In  this  section  the  model  investigated  by  Chomaz  et  al.  (1988)  is 
reexamined  in  the  context  of  the  present  rational  asymptotic  analysis  starting 
from  the  governing  equation  (2.6).  The  main  assumption,  besides  the  exclusion 
of  long-range  feedback,  is  that  the  flow  is  most  (absolutely)  unstable  at  the 
boundary,  i.e.  that  the  location  of  the  "wave  maker"  is  given  by  (3.14). 
Furthermore  we  will  assume  that  3jjU^(X*’)f^0  ,  i.e.  that  Wp(X)  has  no  saddle 
point  at  the  boundary  X‘’-0,  and  that  the  global  mode  amplitude  is  zero  at  the 
boundary . 

The  analysis  is  very  similar  to  the  one  of  section  1.4.  and  will  be  kept 
as  brief  as  possible.  In  the  present  case  the  absolute  growth  rate  <*)^(X)  is 
assumed  to  be  a  linear  function  of  (X-X^)  and  X*'  is  a  first  order  turning 
point.  Hence,  the  balance  between  3^*  and  (X-X*')A*  near  X*'  leads  to 
the  scaling 

X  »  .  (5.1) 

The  disturbance  stream  function  is  expanded  accordingly 
9'  -  ($0  +  +  £2/3$^  +  +  0(£*'®)](X,y)  X 

exp[i£'V(X-X‘')  -  iwgtj  , 


(5.2) 
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where  and  the  global  frequency  is  represented  by 

Wjj  -  ^  Wj  +  euj  +  0(€*^®)  .  (5.3) 

Analugous  to  (4.4),  the  mean  flow  is  expanded  around 
U(X,y)  -  U(X^y)  +  ,y)]  +  0(£*/®)  . 

(5.4) 

V(X.y)  ~  V(X^y)  +  . 

This,  together  with  the  transformation  of  x-derivatives  according  to 

yields,  at  leading  order,  the  same  Rayleigh  problem  (4.5)  at  X^ 
and  the  leading  order  inner  solution  is  thus  determined  up  to  a  free  amplitude 

\(X) 

%(X,y)  -  \a)  4>\i.y)  .  (5.5) 

At  the  next  order  0(£^^®)  one  again  obtains 

-  0  .  (5.6) 

At  order  0(«^^®)  one  finds  the  solvability  condition 

a|So  +  S,{S2  -  XwJ)  -  0  (5.7) 

X^(X— e’^'^X*^)  -  A^(X-^)  -  0  . 

The  coefficient  is  thereby  given  by  (4.10a)  and 

.  (5.8) 

The  first  boundary  condition  on  the  X.  in  (5.12)  expresses  the  assumption 
that  the  global-mode  amplitude  is  zero  at  the  upstream  flow  boundary  X-0, 
while  the  second  boundary  condition  requires  the  amplitude  to  vanish  far 
downstream  as  before.  Equation  (5.7)  is  recognized  as  Airy's  equation  with  the 
solution 

Xj,  -  Ai{(2«5/«^)''MX  -  (V“x>)J  •  (5.9) 

The  boundary  condition  at  X-0  then  leads  to  the  relation  between  and  X*" 

-(2a>$/a)^)^^®  (S/wJ)]  -  -a„  ,  (5.10) 

where  the  -a  are  the  zeros  of  the  Airy  function  Ai  (ap-2.338,  aj-4.088, 
etc.).  With  t^e  location  of  the  tuniing  point  given,  to  leading  order,  by 
(3.14),  i.e.  coinciding  with  the  boundary,  we  obtain 

X*^  -  0  ,  (5.11a) 

The  quantization  of  the  global  frequency  is  therefore  of  order  ,  larger 
than  the  0(e)  quantization  (4.15a)  in  the  doubly- infinite  case.  For  the 
matching  to  the  WKBJ  solution,  it  is  more  convenient  to  have  w^-O.  This 
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can  be  achieved  by  moving  the  turning  point,  aroxind  which  all  qxiantities  are 
expanded,  slightly  away  from  the  boundary.  Setting  S^-O  in  (5.10)  yields 

“2-0  ,  (5.12a) 

K  “  \  •  (5.12b) 

With  the  assumptions  of  a  large -wave -number  cutoff  [i.e.  Im(w^)<0]  and 
maximum  absolute  growth  rate  at  the  boundary  [i.e.  Im(«^)<0]  equation 
(5.12b)  places  the  "dominant  oscillator"  close  to  the  origin  into  the  first 
quadrant  of  the  complex  X-plane.  This  displacement  reflects  the  fact,  that  a 
finite  region  of  absolute  instability  is  required  on  the  real  X-axis  near  the 
origin  before  one  can  have  global  instability,  as  discussed  in  detail  by 
Chomaz  et  al.  (1988)  (see  also  figure  Ic). 


In  the  Fourier -domain  we  again  obtain  the  Taylor  expansion  of  the 
dispersion  relation  around 

€  «3  -  -  e  iw*"  +  (w^/2)(k-k^)2  +  tOx(X-X^) 

(5.13) 

+  (o>^/6)(k-k‘)®  +  a)^(k-k^)(X-X^)  , 

where  it  is  understood  that,  for  the  n-th  global  mode,  all  constants  are 
evaluated  at  X*’-X*’  ,  and  where 

n 

+  3l4(4)  -  61^(4)  -  L^(^^)](L^(^^))-'  . 

(5.14) 

In  these  expressions  the  results  up  to  0(e)  have  been  incorporated  which  is 
carried  out  in  the  preprint  listed  in  section  1.1. 


In  the  two  generic  cases  analyzed  in  this  paper,  in  which  the  global 
instability  is  dominated  either  by  a  saddle  point  of  the  absolute  frequency 
w^(X)  within  the  flow  or  by  one  streamwise  boundary  of  the  flow  domain,  we 
have  developed  approximate  expressions  for  the  global-mode  frequency  its 
growth  rate  and  its  streamwise  amplitude  distribution  in  terms  of  local 
stability  properties  alone.  The  complex  frequency  in  particular  can  be 
estimated  very  easily  from  the  knowledge  of  the  local  absolute  frequency  w 
and  absolute  wave  number  k  on  the  real  x-axis,  except  for  the  nonparallel 
frequency  shift  6u  (equs.  4.15a  and  5.13)  of  order  0(t).  All  that  is  required 
is  the  analytic  continuation  of  w^(X) ,  and  t<>j.(X:kjj)  from  the  real  x- 

axis,  on  which  the  mean-flow  data  are  generally  d^ined,  to  the  respective 
turning  points  X*'  (equs.  3.13  and  5.12b).  Moreover  we  can  conclude  that,  in 
the  context  of  the  present  analysis,  global  modes  only  become  amplified  if  the 
streamwise  extent  of  absolute  instability  is  sufficiently  large:  in  the 
doubly- infinite  case  the  interval  AX  of  absolute  instability  has  to  be  at 
least  of  order  when  X‘'  is  within  of  the  real  axis,  and  larger  for 

X*’  far  from  the  real  axis  (see  Hunt  &  Crighton,  1991).  In  the  semi-infinite 
case,  on  the  other  hand,  global  instability  results  whenever  the  interval  of 
absolute  instability  at  the  end  grows  to  AX-0(€^^^)  (see  equ.  5.12b). 


The  question  naturally  arises  of  how  relevant  our  two  generic  cases  are 
to  practical  applications.  It  can  be  answered  quite  convincingly  for  the  wake 
behind  a  rectangular  plate.  Hannemann  &  Oertel  (1989)  provide  in  their  figure 
16  a  comparison  between  the  linear  from  their  fully  nonparallel  numerical 
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simulation  and  the  local  absolute  frequency  ,  which  clearly  displays  a 

saddle  at  x«l  on  or  very  close  to  the  real  x-axis.  Despite  the  proximity  of 
the  trailing  edge,  it  appears  from  this  comparison  that  the  linear  global 
mode,  here  the  nascent  K4rm4n  vortex  street,  is  dominated  by  a  "wave 
maker"  at  x«l  !  Hence  the  doubly- infinite  model  of  section  1,^.  will  be  used 
in  section  4  to  interpret  the  wake  control  experiments. 


1.6.  Appendices 

A.  The  Rayleigh  operator  and  its  derivatives 

In  the  following  the  Rayleigh  operator  r[»;k,w,X]  and  its  formal 
derivatives  with  respect  to  the  parameters  k,  u  and  X  are  compiled.  In 
addition,  the  operator  which  contains  the  nonparallel  teims  of  the 
disturbance  equation  (2.6)  is  listed  as  (A. 8).  In  all  the  argument  lists  the 
relevant  parameters  are  given  after  a  semi-colon.  For  the  operators  they 
include  the  eigenvalue  pair  k  and  u>  associated  with  a  parallel  flow  U(y;X) 
that  coincides  with  the  local  velocity  profile  'J(X,y)  at  X. 


^(•;k,w,x) 

M 

[kU(y:X)  -  0>]ial  -  k^].  -  k(35u(y;X)]. 

(A.l) 

Xk(*:k.w.X) 

m 

2ka)  .  +U(y:X)(35  -  -  [3jU(y;X)]. 

(A. 2) 

i„(*:k,w,x) 

m 

-[aj  -  k^]. 

(A. 3) 

/jj(*  :k,«,x) 

m 

k(3jjU(y:X)](32  -  k^J.  -  k(3jj3jU(y:X)]. 

(A. 4) 

m 

2u  •  -  6kU(y:X)» 

(A. 5) 

jTjjj  ( • ,  k ,  (a)  ,  X ) 

m 

(ajjU(y:X)]l35  -  3k2).  -  [3^35u(y:X)]. 

(A. 6) 

X^(*  ,k,co,X) 

m 

k(32u(y:X)](32  -  k^).  -  k(3j3y2u(y:X)). 

(A. 7) 

*^vvv  ( *  1  k ,  w ,  X) 

m 

-6U(y:X). 

(A. 8) 

X,(-;k,w.X) 

m 

V(y;X)(35  -  V.\]-  +  [3jj3yU(y  ;X)  )3y.  -  R-'(3j 

-  k^]^*  (A. 9) 

B.  The  k“derivative  of  the  Rayleigh  eigenfunction 

The  derivative  of  the  Rayleigh  equation  (B.l)  with  respect  to  the 
wavenumber  k,  where  £  is  given  by  (A.l),  leads  to  the  inhomogeneous  Rayleigh 
equation  (B.2). 


£i4>^‘,k,o>,X)  -  0  ;  |0o(|y|-^;X)|  -  0  (B.l) 

/(3j^^o:k,w,X)  -  -/j^(^p:k,w.X)  -  /^(^p;k.a),X)  (B.2) 


The  solvability  of  (B.2)  requires  that  the  right-hand  side  be  orthogonal  to 
the  homogeneous  solution  of  the  adjoint  Rayleigh  equation,  i.e.  to 
[kU(y:X) -w]‘^  With  the  notation  (B.3)  for  the  inner  product 
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«|  jr  («;k,a),X)  -  dy  ,  (B.3) 

yi  kU(y;X)-« 

onj?  obtains  a  convenient  expression  (B.4)  for  the  complex  group  velocity 
in  terms  of  the  Rayleigh  eigenfunction 

[L„(^o:k,w,X)]‘^  -  0  (B.4) 


C.  The  X-derivative  of  the  Rayleigh  eigenfunction 


In  an  analogous  manner  an  expression  for  the  X-derivative  of  cau  be 
obtained.  Differentiation  of  (B.l)  with  respect  to  X  leads  to  the 
inhomogeneous  Rayleigh  equation  (C.l>. 


£id^<f>^]k.w,X)  -  -a^jk  JCj^(^Q:k,w.X)  -  8^  -  /jj(^g:k,«.X) 

With  the  notation  (B.3),  the  solvability  condition  reads 

ajjk  i^(^Q:k.w,x)  +  a^Lj^^‘,k.u>,x)  +  i^(^^;k,u),x)  -  o  . 

Using  the  result  (B.4),  this  can  be  recast  in  the  form 
-  3jjk  d^u>  +  a^  +  l^(^^-,k,o>,X)  (L^(^g;k,6),X)]''  -  0  . 

For  use  in  §3  we  also  list  a  particular  solution  of  (C.l) 


(C.l) 


(C.2) 


(C.3) 


Vo  “  'V  ^uc  V  K  ■ 

in  which  the  functions  4>^(y',X) ,  etc.  are  particular  solutions  of  the 
inhomogeneous  Rayleigh  equations  (C.5)  which  satisfy  the  same  boundary 
conditions  as 


/(^^j:k,a),x)  -  i:j(^g:k,«,x)  ;  !^i,(ly|^:X) |  -o  .  (c.5) 

In  addition  we  will  also  need  the  following  "second-generation"  forced 
solution 


£i<i>^-,k,u>,x)  -  j:^(^jj^;k,w,x)  :  i^2j^(|y!->«>:X)|  -  o 

where  is  given  by  (C.5). 


(C.6) 
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2.  SINGLE-INPUT  SINGLE-OUTPUT  CONTROL  OF  GLOBAL  MODES;  REDUCTION 
OF  THE  PLANT  MODEL  TO  A  STUART-LANDAU  (SL)  EQUATION,  THE 
CONTROLLER  DESIGN  AND  ITS  IMPLEMENTATION  IN  A  HEATED  2-D  JET 

2.1.  Related  Publications  and  Presentations 

HENRI CH,  EDWARD  A. ,  Control  of  limit  cycle  oscillations  with  applications  in 
fluid  mechanics.  Fh.D.  Thesis,  University  of  California,  Los  Angeles,  1991. 

HENRICH,  EDWARD  A.,  MINGORI,  TINO  D.L.  &  MONKEWITZ,  PETER  A.,  Control  of  Hopf 
bifurcations.  Proc.  of  NATO  ASI  on  Chaotic  Dynamics:  Theory  and  Practice, 
Patras,  Greece,  July  1991. 

GAMET,  LIONEL,  Control  of  a  two-dimensional  hot  Jet.  DEA  Thesis  (ENSICA, 
Toulouse)  carried  out  at  the  University  of  California,  Los  Angeles,  1991. 

MONKEWITZ,  PETER  A.,  GAMET,  LIONEL,  HENRICH.  EDWARD  A.  &  MINGORI,  TINO  D.L. , 
Feedback  control  of  a  self-excited  heated  two-dimensional  jet.  Bull  Am.  Phys. 
Soc.  Vol.  36.  p.  2629,  1991. 

MONKEWITZ,  PETER  A.,  MINGORI,  TINO  D.L. ,  HENRICH,  EDWARD  A.  &  YU,  MING-HUEI, 
Adaptive  and  nonadaptive  feedback  control  of  global  instabilities.  Proc.  AFOSR 
Contractors  Meeting  on  Turbulence  -  Structure  and  Control,  Columbus  Ohio,  pp. 
103-106,  April  1991. 

HENRICH,  EDWARD  A.,  MINGORI,  TINO  D.L.  &  MONKEWITZ,  PETER  A.,  Control  of 
Pitchfork  and  Hopf  bifurcations.  Submitted  to  1992  ACC. 

HENRICH,  EDWARD  A.,  MINGORI,  TINO  D.L.  &  MONKEWITZ,  PETER  A.,  Control  of  limit 
cycle  oscillations  with  applications  in  fluid  mechanics.  To  be  submitted  to 
Int.  J.  Control,  1992. 


2.2.  Introduction  and  the  Analysis  of  a  Hopf  Bifurcation  by 
Multiple  Scales 


Many  systems  in  nature  develop  self-excited  or  limit  cycle  behavior  when  the  most  amplified  global 
mode  becomes  temporally  unstable.  In  fluid  mechanics,  two  examples  are  a  wake  behind  a  cylinder 
and  a  heated  jet  [1,2].  The  goal  of  this  research  is  to  control  the  self-excited  oscillations  exhibited 
by  such  systems  or  to  induce  them  in  globally  stable  systems. 

For  the  control  of  flow  oscillations,  the  primary  task  is  the  development  of  simplified  model  as  the 
underlying  partial  differential  equations  are  too  unwieldy.  In  Section  1  we  have  derived  lineariied 
GL  equations  describing  global  modes  in  doubly-  aad  semi-infinite  systems  which  can  be  extended 
to  the  weakly  nonlinear  regime.  When  considering  only  a  single  gobla  mode,  the  GL  equation 
can  be  reduced  to  a  Stuart-Landau  equation  for  its  characteristic  amplitude,  which  is  an  ordinary 
differential  equation  in  time.  The  latter  can  be  viewed  as  describing  the  output  of  the  probe  at 
some  fixed  space  location.  The  following  development  is  carried  out  in  this  spirit  and  the  resulting 
mathematical  model  for  control,  not  being  specific  to  the  wake  or  the  hot  jet,  should  be  abie  to 
control  limit  cycles  in  a  variety  of  physical  systems. 

A  common  bifurcation  in  fluid  naechanics  is  the  Hopf  bifurcation  of  a  steady  state  to  a  limit  cycle. 
In  many  cases,  one  bifurcates  from  the  tero  state  to  a  non-iero  steady-state  to  a  limit  cycle  as  a 
parameter  is  mcreased  These  bifurcations  are  generally  the  first  steps  in  the  transition  to  chaos  or 
turbulence  [3,  4]  For  a  given  physical  system  there  are  often  design  constraints  which  specify  a  fixed 
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region  of  pariumter  space  in  which  one  would  like  to  operate.  If  this  particular  chcnce  oi  parameters 
exhibits  some  undesirable  dynamic  behavior  (e.g.  limit  cycles  or  chaos)  then  the  goal  of  the  control 
engineer  is  to  suppress  this  behavior  while  still  operating  in  the  same  parameter  range.  A  distinction 
is  made  between  a  ‘control  parameter’  or  ‘bifurcation  parameter’  which  is  the  parameter  <Hie  uses 
to  study  successive  bifurcations  and  a  ‘control  force’  which  is  a  prescribed  function  of  time  and 
possibly  space  used  to  modify  the  system’s  output.  As  an  example,  consider  that  one  might  desire 
to  suppress  the  Von  Karman  vortex  street  in  the  wake  behind  a  cylinder.  The  control  parameter  is 
the  Reynolds  number  which  may  be  fixed  at  some  supercritical  value,  where  the  wake  is  no  longer 
steady.  The  control  force  could  be  an  induced  acoustic  field  designed  to  eliminate  flow  oscillations. 
This  method  uses  the  control  force  to  stabilize  the  system  in  contrast  to  methods  which  modify  the 
bifurcation  parameter  until  the  state  behaves  as  desired  [5]. 

In  an  attempt  to  understand  the  underlying  properties  of  control,  the  research  reported  here 
concentrates  on  the  development  of  control  laws  for  finite  dimensional  nonlinear  systems  which 
undergo  a  flopf  bifurcation  to  a  limit  cycle.  The  relevance  to  fluid  oscillations  is  that  the  finite 
dimensional  nonlinear  system  can  be  thought  of  as  a  finite  difference  or  Galerkin  approximation  to 
the  Navier-Stokes  equations.  The  Stuart-Landau  equation  is  valid  when  the  eigenvalues  of  these 
approximations  are  well  separated. 

Usmg  a  multiple  scale  perturbation  technique  we  derive  a  Stuartr Landau  equation  which  governs 
the  effect  of  a  linear  regulator  on  a  system  undergoing  a  bifurcation.  Often  one  does  not  know  the 
differential  equations  governing  the  system,  but  one  does  know  that  the  system  exhibits  behavior 
indicative  of  a  Hopf  bifurcation.  In  this  case  one  still  knows  that  asymptotically  the  dynamics  are 
governed  by  a  Stuart-Landau  equation,  so  that  this  simplified  equation  can  be  used  as  a  generic 
mode!.  One  can  then  estimate  the  parameters  of  the  Stuart-Landau  equation  using  experimental 
data  and  then  determine  an  appropriate  control  strategy  based  on  this  information. 

The  Hopf  bifurcation  is  a  branching  of  time-periodic  solutions  firom  an  equilibrium  solution  branch. 
When  the  bifurcation  is  supercritical,  as  one  increases  a  parameter  /i  the  ^uilibrium  solution  be¬ 
comes  unstable  and  a  limit  cycle  is  born.  As  one  increases  the  control  parameter  further  one  secs  more 
and  more  complicated  behavior.  These  successive  bifurcations  comprise  the  basis  for  several  theories 
for  the  transition  to  turbulence  including  period  doubling,  intermittency  and  quasi-periodicity  (3, 4]. 
We  now  derive  an  amplitude  equation  for  the  Hopf  bifurcation,  the  development  follows  that  of  [6]. 

For  a  physical  system  with  a  finite  number  of  degrees  of  freedom,  the  dynamics  may  be  descried 
by  a  set  of  first  order  ordinary  differential  equations  in  vector  fo  ' 

^  =  F(X(0;p)  (2.1) 

where  X  e  R"  is  the  system  stale,  F  is  a  vector  field  and  /r  €  R  is  a  control  parameter  in  the  problem. 
Equation  (2.1)  is  presumed  to  have  cn  equilibrium  solution  Xo(/r)  which  satisfies  F(Xo(pj;/i)  =  0. 
We  can  then  express  (2.1)  as  a  Ihylor  series  about  the  equilibrium  point  Xo(/i)  so  that 

dx 

—  =  Lx-h  Mxx  -f  Nxxx  -f  •  •  •  (2.2) 

where  the  expansion  is  in  terms  of  x  =  X  -  X©.  For  the  Thy  lor  series  to  be  valid  we  need  x  to  be 
small,  so  that  the  nonlinear  terms  in  (2.2)  are  small  compared  to  the  linear  terms,  i.c.  the  system 
is  “Sveakly  nonlinear”.  The  matrix  L  is  the  Jacobirm  whose  elements  are  given  by 

Q 

=  ^^(^;)  U=Xo  (2.3) 

The  abbreviations  Mxx  and  Nxxx  denote  vectors  coistaining  the  quadratic  and  cubic  terms  in  the 
expansion.  The  jth  component  of  these  vectors  is  given  by 

(Mxx),  =  ^  X! 

k  I  fn 


(2.5) 
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We  are  interested  in  Hopf  bifurcetions,  so  we  introduce  some  restrictions  on  the  linear  operator 
L.  The  stability  of  the  eqidlibrium  solution  Xo  is  determined  by  the  eigenvalues  of  L.  We  msuttw 
that  L  is  stable  for  p  <  0  and  unstable  for  ft>Q.  By  L  unstable,  we  mean  that  at  least  one  of  the 
eigenvalues,  A,  fcimd  by  solving  Lr  =  Ar  has  a  real  part  greater  than  zero.  We  also  M«uTn«»  that 
these  eigenvalues  have  nonzero  transversal  ‘Velocity”  when  crossing  the  imaginary  axis.  That  is 


dfi 


b.=:O>0 


(2.6) 


which  we  will  shortly  see  fixes  the  scaling  for  the  slow  time  scale  to  be  |/i  1 1. 
We  expand  the  linear  operator  in  powers  of  fi 


(2.7) 


sinoilar  expressions  may  be  developed  for  M  and  N.  The  eigenvalues  of  X  are  expressed  as 

Ai(p)  =  Af +  /«A<‘>  +  /,’Af>+...  (2.8) 

The  convention  that  a  superscript  in  parentheses  denotes  a  perturbation  and  a  subscript  refers  to 
an  index  will  be  used  throughout.  The  eigenvalues  are  complex  m  general  and  can  be  written  as 
A^*^  =  Let  tj  and  1^  denote  the  right  and  left  eigenvectors  of  X^®^  so  that 


=  xf^Tj  (2.9) 

IjXW  =  Xfhj  (2.10) 


For  simple  eigenvalues,  the  eigenvectors  are  orthogonal  (l^rt  =  0  for  j  ^  k)  and  we  can  chooae  to 
normalize  such  that  l^r^'  =  1.  We  now  define  a  small  positive  pvanoeter  <  such  that  fi  ~  where 
X  s=  8gn(/i).  We  can  expand  the  state  in  powers  of  «  as 

X  =  +  •  •  •  (2.11  y 

Note  that  the  expansion  starts  with  terms  of  order  e.  This  follows  from  the  assumption  tha>  (he 
system  (2.1)  is  weakly  nonlinear.  For  systems  where  this  is  not  the  case  the  equations  {2.’’ 6)  are 
not  generally  solvable  since  the  0(1)  equation  is  also  nonlinear.  Also  note  that  defining  t  —  \/\ti\ 
represents  a  balance  of  the  relevant  nonlinear  terms  and  the  linear  dampmg  or  dedampmg  of  the 
critical  modes. 

Since  the  real  part  of  the  critical  eigenvalue  is  O(f^)  it  is  natural  to  mtroduce  a  second  time 
scale  T  =|p  1 1  =  We  treat  the  two  time  scales  as  independent  variables  so  that  the  x^^l’s  are 
functions  of  two  variables  [7].  The  rule  for  calculating  derivatives  becomes 


di 


L 

dr 


(2.12) 


Inserting  these  expansions  in  p  and  e  into  the  Tzylor  expansion  (2.2)  and  equating  like  powers  of  < 
results  in  an  infinite  set  of  linear  inhomogeneous  partial  differential  equations  of  the  form 


(”-X<°))xW  =  ^‘^  (2.13) 

The  first  few  are 

0(e)  :  =  0  (2.H) 

O(e^) :  =  MWx(‘)x(i)  (2.15) 

0(c*) :  =  2M(0)x(‘)x(*>  +  iVWx<‘)x(5)x<‘>  -  (^  -  A’X<‘))x<‘>  (2.16) 


Note  that  the  are  functions  of  x^^  with  j  <  i  bo  that  they  are  known  functions  if  we  solve 
the  set  of  equations  in  order. 
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The  lowest  order  problem,  0(e),  is  linear  and  homogeneous,  and  the  higher  order  problems  are 
linear  and  inhomogeneous.  The  0(e)  equation  is 

(2.17) 

at 

which  has  the  solution  in  terms  of  the  eigenvectors 

^^%.r)=^Ci(T)t^f\  (2.18) 

i=i 

Since  (2.17)  is  a  partial  differential  equation,  the  solution  is  found  in  terms  of  an  arbitrary  function 
of  the  slow  time  scale  r.  For  this  reason,  the  coefficients,  Cj ,  are  aUowed  to  depend  on  the  alow  time 
scale  T  —  e^t  in  standard  multiple  scale  fashion.  The  hi^er  order  problems  determine  conditions 
on  these  arbitrary  functions  in  order  to  make  the  solution  uniformly  valid  for  times  of 

0(l/e).  One  must  impose  “solvability”  conditions  so  that  the  solution  x^^^(f,T)  does  not  involve 
any  “secular”  terms  (i.e.  terms  which  become  unbounded).  Since  the  operator  on  the  left  hand  side 
of  (2.13)  is  the  same  for  each  order,  the  solvability  condition  is  that  we  must  eliminate  terms  on  the 
right  hand  side  of  (2.13)  which  are  marginally  stable  solutions  of  the  hunogeneous  equation. 

We  number  the  two  critical  eigenvalues  first  so  that  and  rewrite  (2.18)  as 

x('>(f ,  t)  =  A(r)e^^%t  +  i(r)c-‘““*‘ri  +  (2.19) 

where  A(t)  =  ci(r)  =  C2(t)  will  be  the  dependent  variable  in  a  complex  amplitude  equation.  Since 
the  real  part  of  the  critical  eigenvalues  is  0(«®)  the  slow  temporal  growth  or  decay  is  absorbed  in 
the  slowly  varying  amplitude  A(t). 

The  second  equation,  O(e^),  is 

^,(3)  ^  +  M(®)x(‘)x(‘)  (2.20) 

To  solve  the  O(e^)  equation  we  note  that  since  is  quadratic  it  acts  like  an  input  forcing  at 
and  D.C.  as  well  as  other  nonsecular  inputs.  The  solution  is 

x<’)  =  V+ylV‘“'‘’^‘  +  V_iV^‘'”‘  +  Vo  |yll’+---  (2.21) 

where  +  •  •  •  are  terms  which  will  not  cause  a  secularity  at  O(e^).  Substitution  of  (2.19)  and  (2.21) 
into  (2.20)  and  equating  coefficients  of  A^,  A^  and  |j41’  and  solving  for  V+,  V_  and  Vo  we  <^tain 

V+  =  V_  =  -  2iw(°)]-*AfWriri  (2.22) 

Vo  =  -2{L<®)]-»M(®)rir,  (2.23) 

Now  we  concentrate  on  elinunating  the  secular  terms  at  order  t*.  lb  do  this  form  “*1 

then  set  the  coefficients  of  e^^***’'  equal  to  sero,  which  gives  the  solvability  condition  at  order 

li{2i^’M<®>rjV+  +  2  li4 1®  ^Af<®)riVo 

+3fV(°)rjrif -  ^Arj  +  XL<^^'>Ati]  =  0  (2.24) 

This  can  be  rewritten  as  an  amplitude  equation 

^  =  X>S^'>A-g\A\^A  (2.25) 


where 


(2^6) 

(2.27) 


g  =  -21iA/(°>fiV+  -  21iAf(°)riVo  -  31i 

Equ&tion  (2.25)  is  knowo  as  the  Stuart-Landau  equation. 

Example:  Kan  der  Pol  Oscillator 

To  illustrate  the  method,  we  consider  the  Van  der  Pol  equation  in  dimensionless  form 

V  +  (»^-/i)y  +  »  =  0  (2.28) 

If  we  take  the  state  to  be  X  =  [y  then  the  vector  field  is 

F(X;rt=[_;  i]x+[.xU] 

The  equilibrium  state  is  Xo  =  [0  0]^,  so  that  the  Jacobian  is 

^=[-1  o]+^[o  +  (2-30) 

and  the  cubic  term  in  the  Thylor  expansion  is 

1 

(JVxyz)2  =  -  j(*2yi2i  +  xig^zi  +  iivizj)  (2.31) 

The  left  and  right  eigenvectors  for  the  critical  eigenvalue  =  i  are  Ij  =  ^(1  —  i]  and  ri  =  [1 

The  perturbation  in  the  eigenvalue  is  Aj'^  =  liL^'^rj  =  The  coefBcient  of  the  nonlinear  term 

in  the  Stuart-Landau  equation  is  y  =  — 31i^(°)ririfi  =  The  Stuart-Landau  equation  for  this 

example  is 

^^IxA-\a\A\^  (2.32) 

The  initial  condition  for  the  Stuart-Landau  equation  is  j4(0)  -  ^(y(0)  -  iy(0)).  Figure  1  compares 
the  exact  solution  of  the  Van  der  Pol  equation  (by  numerical  integration)  to  the  asymptotic  solution 
foimd  by  solving  the  Stuart-Landau  equation.  In  this  figure  our  small  parameter  is  (  =  .5.  llie 
important  thing  to  note  is  that  the  short  term  prediction  using  this  model  is  very  good,  and  short 
term  prediction  is  all  that  is  needed  for  control,  n 

2.3.  The  Linear  Regulator  and  the  Nonlineeir  Plant 

This  paper  addresses  the  problem  of  controlling  a  system  by  external  forcing  when  the  control 
parameter  b  fixed  by  some  other  design  constraint.  We  assume  that  the  plant  (subscript  p)  b  a  single¬ 
input  single-output  system  with  nonlinear  dynamic  equations  which  undergo  a  Eopf  bifurcation  at 

p  =  0 


^  =  F(X,:p)-fB,u 

(3.1) 

y  =  C,X,  +  D, 

(32) 

X,6iZ"'  and  /:,u€R 

(32) 

The  goal  of  the  control  system  h  to  return  the  system  to  its  equilibrium  state  which  has  become 
jpen  irop  unstable  for  ft  >  0,  A  nrAural  choice  for  a  control  system  b  the  linear  regulator  as 
depicted  in  Figure  1.  The  differcutial  equations  governing  the  control  syi  iem  are  such  that  the 
Laplai.e  transform  of  impulse  .■-sponse  of  the  linear  compensator  (subscript  c)  b  KG{$). 

■'he  eoja.ions  govermag  the  regulator  are 


Figure  1:  Ckimparison  of  Exact  and  Asymptotic  solutions  of  the  Van  der  Pol  equation  for  e  =  .5  a) 
Phase  plane  for  exact  solution  b)  Phase  plane  for  asymptotic  solution  c)  Exact  solution  d)  Difference 
between  exact  and  asymptotic  solutions 


Figure  2:  Block  diagram  of  the  control  system 
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Xe  =  AeXc  +  KBgt 

(3.4) 

tl  =  CeXc  +  KDet 

(3.5) 

€  and  e,/i  6  R 

(3.6) 

Where  e  =  r  —  y  is  the  error  signal,  and  the  reference  signal  is  the  output  which  corresponds  to  the 
equilibrium  solution  r  =  yo  =  CpXpo  +  Dp.  The  compensator  is  asymptotically  stable,  that  is,  the 
eigenvalues  of  At  lie  in  the  open  left  half  plane. 

We  may  expand  (3.1)  in  a  Taylor  series  about  its  equilibrium  point  X,k)  so  that 

=  LpXp  +  MpTLpXp  +  NpXpXpXp  f  — b  Bpti  (3.7) 


where  F(XpOi  /<)  =  0  defines  the  equilibrium  point  and  Xp  ~Xf  —  X^o  is  the  expansion  variable. 
The  output  equation  becomes  y  =  CpXp  +  CpXpo  +  Dp.  Combining  these  two  vector  equations  bto 
a  composite  equation  we  obtidn 


^  r  Xp 
dt  Xc 


Lp  -  KBpDcCp  BpCc  1  [  *!>  1  a.  f  +  NpXpXpXp  +  •  •  • 

~KB,Cp  J  I  Xc  J  [  0 


which  is  in  the  standard  form  (2.2).  We  assume  that  the  plant  undergoes  a  Eopf  bifurcation  at 
/i  =  0  and  that  we  can  expand 

ip  =  4°>  +  /i4*)  +  .--  (3.9) 


in  the  usual  way.  We  assume  that  we  choose  a  small  gain  such  that  KCp  ~  0(/i)  and  so  we  rewrite 
K  =1/^1  K.  The  expansion  for  the  composite  system’s  linear  operator  is 


where 


and 


i(o)  _  [  BpCc 

0  Ac 


(3.10) 

(3.11) 


(3.12) 


£,(1)  _  -  Bga{ti)KBpDcCp  0 

-8gn(/i)/cBeCp  0 

Order  the  plant’s  eigenvalues  such  that  the  marginal  eigenvalues  are  first  and  denote  the  plant’s 
left  and  right  eigenvectors  by  lipLp°^  =  ta;(°lljp  and  Lp’^rip  =  The  right  eigenvector  for 

the  composite  naatrix  is  ri  =  [  rfp  Ow,  and  the  left  eigenvector  is 

li  =  (  lip  lipBpC.(twW-Ac]-»  ]  (3.13) 

We  compute  the  perturbation  in  the  critical  eigenvalue  as 

=  liL(‘)ri  =  lip4‘)r,p  -  sgn(^)lipBC7rip»c[Cc(ia-W  -  ^c)~‘Bc  +  £>c)  (3.14) 


We  can  rewrite  (3.14)  as 

Ai‘^  =  A(‘>-A'KJipG(tw(‘»)  (3.15) 

where  is  the  uncontrolled  perturbation  in  the  critical  eigenvalue,  Jip  =  lipBCrip  is  the  recep¬ 
tivity  of  the  critical  mode  and 


KG{8)  =  K{Cc{tI  -  Ac)-^Bc  -f  Dc] 


(3.16) 
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is  the  tr&nsfer  function  of  the  linear  compensator.  We  see  that,  in  the  context  of  multiple  scales,  all 
compensators  with  the  same  gtun  and  phase  at  the  critical  fi^uency  have  the  same  effect  on  the 
plant.  Note  that  Jip  is  what  is  known  as  a  "modal  inSuence  coefficient”  in  structural  mechanics 
or  the  "receptivity”  in  fluid  mechanics  corresponding  to'the  critical  mode.  This  complex  constant 
determines  how  the  control  input  couples  to  the  output.  The  goal  of  the  control  system  is  to  force 
"in  phase”  to  enhance  oscillation  and  "out  of  phase”  to  reduce  oscillation.  The  coefficient  Jxp 
us  which  phase  is  "in”  and  which  phase  is  "out”. 

The  nonlinear  term  in  the  Stuart-Landau  equation  is  unchanged  by  the  addition  of  this  linear 
compensator.  The  Landau  constant  remains  g  =  gp.  To  see  this,  first  compute  0;vc 

and  Vj  =  [  V^,  Oatc  ]^*  Then 

g  =  -21ipMWripV+,  -  aipM^npYop  -  31,pN;'’)ripr,pfip  =  Pp  (3.17) 

Hence,  the  Stuart-Landau  equation  for  the  regulated  system  is 

^  -  JipKG{iu,(%A-gp  A  (3.18) 

The  preser\'ation  of  the  nonlinear  term  is  a  consequence  of  choosing  the  compensator  gain  to  be 

0(/i). 

Example:  Visn  der  Pol  Oscillator  toilh  a  Linear  Regulator 
The  equation  for  the  Van  der  Pol  oscillator  with  a  control  force  is 

y+(y^-p)v+v  =  «  (3.19) 


The  StuartrLandau  equation  for  the  regulated  system  is 


^  =  l(X  +  iKC{i))A-hAfA 


(3.20) 


The  method  of  multiple  scales  predicts  the  closed  loop  system  will  be  stable  fot Im(KG{i))  >  l.If 
we  consider  derivative  feedback  so  that 


u(«)  =  KG{s)v{s)  =  Ksv[s) 


(3.21) 


then  we  expect  stability  for  if  >  ^.  It  is  easy  to  check  that  this  is  in  fact  the  exact  result.  If  we 
consider  a  lead  compensator  so  that 

=  (3.22) 

then  the  multiple  scales  method  predicts  closed  loop  stability  for 

K{b  —  a) 


For  p  >  0  the  closed  loop  system  is  stable  for 

it(6  —  a) 


>  1 


62 -hi 


>  1-f  p 


K-b 
62 -fl 


(3.23) 


(3.24) 


So  the  multiple  scales  stability  result  is  correct  to  0(1). 

For  comparison,  we  consider  the  same  plant  and  control  and  determine  the  center  manifold 
reduction  and  normal  form  [8].  The  dynamics  are 


p  =  0  (3.25) 

V  =  v  (3.26) 

v  = -y  +  pv  — y’v-bu  (327) 

XI  =  ~b\i  -  kX fiv  —  kX ang  (328) 
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where  terms  involving  fi  are  nonlinear  since  ft  is  also  a  state  variable.  We  have  a  three  dimensional 
center  manifold.  The  center  manifold  is 


,  .  .„l  +  a5  ■  ^b  — a 

u  =  u{ft,  y,  v)  =  -kX  jjqjTj  PV  -  kX  +  •  •  • 

The  dynamics  on  the  center  manifold  are  governed  by  the  equation 

r«^(h-a)  ,1.  r  ic;i'(l  +  a6)1 

"+»=-'■  l-jf+i  ■  ^ '  r  TO--^J ■ 

Reducing  this  equation  to  normal  form  gives 


#c£(6-^1  1, 

dt  2[  fc’  +  l  J“8 


(3.29) 


(3.30) 


(3.31) 


which  is  identical  to  the  multiple  scales  result  (^4  =  re**/(2c)tT  =  ^i).  O 

By  specifying  the  gain  of  the  compensator  to  be  0{fi)  we  have  derived  an  amplitude  equation 
which  explicitly  shows  the  effect  of  any  linear  compensator  on  a  nonlinear  system  near  a  Hopf 
bifurcation.  The  importance  of  this  result  goes  beyond  analyzing  the  effect  of  a  control  system  on 
the  plant.  Since  any  plant  which  exhibits  the  features  of  a  Hopf  bifurcation  is  asymptotic  to  the 
solution  of  some  Stuart-Landau  equation,  we  can  use  the  Stuart-Landau  equation  as  a  generic  model 
for  a  Hopf  bifurcation  even  when  the  plant  equations  are  unknown.  We  can  assume  the  plant  is  the 
Stuart-Landau  equation  and  then  try  to  estimate  its  coefficients.  We  may  select  several  compensators 
in  order  to  estimate  the  coefficients  and  Ji,  through  control  induced  transients.  Once  we  have 
accurate  estimates  for  these  parameters  we  may  specify  the  correct  phase  of  the  compensator  at  the 
aitical  frequency.  That  is  we  can  determine  the  best  G(tu)(°l). 

Using  the  same  type  of  analysis  we  find  that  equation  (3.18)  holds  for  pitchfork  bifurcations  as 
well  if  we  set  =  0  (9, 10]. 


2.4.  An  Amplitude  Equation  for  Pareuaeter  Estimation 

To  control  a  given  physical  system,  we  intend  to  use  a  feedback  control  based  on  a  set  of  measure¬ 
ments  which  depend  on  the  system  state.  A  general  nonlinear  measurement  equation  y  =  H{X)  can 
be  written  as  a  Taylor  series 

V  =  -f  -I-  -!-•••  (4.1) 

Since  x  is  O(t'),  the  leading  order  contribution  is  the  affine  equation 

y=CX  +  D  (4.2) 

The  leading  order  effect  of  the  measurement  equation  is  to  scale  the  system  state  and  provide  a 
D.C.  bias.  If  the  system  we  are  interested  in  contains  a  small  parameter  ft^  X^  and  exhibits  limit 
cycles  when  p  >  0,  we  detect  c  not  by  the  size  of  the  output,  but  by  the  slow  time  scale  and  by  how 
far  the  phase  plane  diagram  of  the  limit  cycle  is  "out-of-round”.  In  the  limit  as  e  0  the  output 
is  a  sifioeoid.  So  that  in  terms  of  the  physical  observation,  the  term  "weakly  nonlinev”  implies  the 
system  Las  “almost  sinusoidal”  oscillations.  The  leading  order  asymptotic  expression  for  x  is 

X  ^  €[A(r)e*“'*’*‘ri  -f  c.c.]  (4.3) 

We  define  Y  (7  )  ai>  iLe  output  “amplitude” 


y(r)  =  2f  A(r)Cri 


(4.4) 
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F^om  equations  (4.2),  (4.3)  and  (4.4)  the  leading  order  asymptotic  expression  for  the  output  is 

V{t)  ~  2iRc{A{T)Ctit^'^^^]  +  CXo  +  2>  =  i2e[y(r)e^‘*’‘]  +  CXo  + 1?  (4.5) 

Equation  (4.4)  relates  the  physical  amplitude  Y (r)  to  the  amplitud'  in  the  normalized  Stuart-Landau 
equation  A{t).  We  are  interested  in  real-time  estimation,  so  we  also  return  to  the  physical  tim<» 
scale  t  =  r/e’.  Making  the  appropriate  substitutions  in  the  Stuart-Landau  equation  (3.18)  leads  to 
an  amplitude  equation  for  y(t) 


where 


■^  =  A^^Y-g\Y\^Y  +  JU(t) 


-21iA/(‘>)fiV+  -  21iM(®)r,Vo  -  31i7V(‘>>r,rir, 

4!CriP 

J  =  (C,r,)(liB,) 


(4.6) 


(4.7) 

(4.8) 

(4.9) 


and  U{i)  =  KG(iu/(^^)Y(t)  for  the  linear  regulator.  Note  that  this  is  an  input-output  description. 
The  coefficient  J  combines  the  efliects  of  control  and  observation  like  a  “loodal  influence  coefficient” 
used  in  model  reduction. 

The  dimensional  output  y(t)  is  recovered  from  (4.5).  The  dimensional  value  of  the  control  is 
given  by 

u«)  =  /Ze(i/(t)e‘-‘“’‘)  (4.10) 

Notice  that  we  have  eliminated  the  unknown  parameter  <. 

The  complex  amplitude  equation  which  models  Eopf  bifurcations  is 

y  =  Ay-p|ypy  (4.ii) 

where  A  =  without  control  and  A  =  =  A°^  -  JKG(vJ-°'>)  with  feedback  control.  One 

can  demodulate  the  measured  output  v(<)  using  a  Hilbert  transform  to  get  Y{i)  (9,  11].  We  consider 
parameter  estimation  for  the  Stuart-Landau  equation  with  transient  data.  We  use  a  Lyapunov 
function  E  =|y  P  to  develop  a  procedure  for  estimating  the  real  parts  of  A®^  and  g.  FVom  (4.11) 
we  have 

E  =  2ArE-2grE^  (4.12) 

If  Ar  >  0  and  yr  >  0  a  stable  limit  cycle  exists.  If  A,.  <  0  and  gr  >  0  then  the  Stuart-Landau 
equation  is  globally  asymptotically  stable.  We  expect  that  a  feedback  which  makes  Ar  <  0  in  the 
closed  loop  should  stabilize  the  system,  that  is  it  should  eliminate  limit  cycle  oscillations. 

Since  we  are  primarily  interested  in  controlling  the  amplitude  of  the  limit  cycle,  it  is  sufficient  for 
our  purposes  to  estimate  the  coefficients  Ar  and  g^  The  coefficients  A^  and  y,  can  also  be  determined 
(see  [11]). 

To  estimate  the  parameters  of  (4.12),  we  can  approximate  the  derivative  in  order  to  estimate  the 
parameters  of  the  Stuart-Landau  equation.  Using  an  Euler  approximation  to  the  time  derivative 
equation  (4.12)  becomes 

E(k  -fl)  =  (1  -f  2ArA)E{k)  -  (2yr  A)£:(k)’  (4.13) 

The  equilibrium  solutions  of  equation  (4.13)  are  E{oo)  =  0  and  E{oo)  =  Ar/yr  which  is  the  same 
as  the  equilibrium  solutions  of  equation  (4.12),  independent  of  the  sample  time,  A. 

So  we  should  be  able  to  replace  the  continuous  time  Stuart-Landau  equation  (4  '  1)  by  an  ap¬ 
proximate  sampled  data  equation  in  the  form 
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E{k+i)  =  eiE{k)  +  e2Eikf 


(4.14) 


Since  this  equation  is  linear  in  the  parameters  we  can  derive  a  least  squares  estimator  for  the 
parameter  vector  0.  As  a  cautionary  note,  we  observe  that  we  have  replaced  the  well  behaved 
continuous  time  logistic  equation  (4.12)  with  a  discrete  logistic  equation  or  quadratic  map  (4.14), 
which  is  known  to  exhibit  chaos  for  a  rimge  of  its  parameters.  For  0<  Ar  <  f$  that  the  dynamics 
Will  asymptote  to  the  fixed  point  E{oo)  =  Ar/9r-  That  is,  as  long  as  the  sampling  rate,  f,  =  1/A, 
is  larger  than  the  slow  exponential  growth  rate,  Ar,  there  will  be  no  periodic  or  chaotic  solutions  to 
equation  (4.13)  [9]. 

F^om  the  batch  form  of  the  least  squares  problem 


(4.15) 


we  could  also  derive  a  formula  for  the  least  squares  parameter  estimate. 

JE7  =|y  P  is  known  since  we  can  find  Y  by  usbg  the  Hilbert  transform.  The  estimation  procedure 
works  best  when  the  initial  conditions  are  small  so  that  we  have  some  data  in  the  linear  growth  range. 
In  practice  this  might  require  several  parameter  excursions  below  the  critical  value  of  /i  to  obtain 
a  steady  state  which  is  a  fixed  point.  Then  we  could  increase  p  above  critical  and  get  a  transient 
with  some  data  in  the  linear  range.  In  principle  the  system  could  be  identified  by  using  parameter 
excursions,  but  for  the  hot  jet  this  is  impractical  since  we  cannot  instantaneously  change  the  density 
of  the  hot  fluid.  A  more  practical  parameter  estimation  scheme  is  to  use  a  feedback  control  system 
to  induce  transients. 

If  we  consider  a  closed  loop  control  with  a  fixed  gain  the  closed  loop  Stuart-Landau  equation  is 

y  =  (A°^  -  jifG(iw<®>))y  -  p  ly  I’  y  (4.16) 

We  can  use  the  same  procedure  on  the  closed  loop  system  as  the  uncontrolled  system  to  find 

=  Af^- 7,^008^  + Jiff  sin  (4.17) 

and  Qr.  Using  “control  induced  transients”  is  a  more  practical  approach  to  estimation.  By  selecting 
several  different  gains  and  phases,  we  can  use  the  data  from  each  of  these  control  bduced  transients 
to  determine  the  parameters  Ar,  Jr  and  J,  by  least  squares. 


1  -/ficoe^r 

/fjsb^i 

\  1 

■  A^;f  ■ 

•  ; 

j 

Jr 

r= 

1  -KflVMiflN 

Jff/vsin^/v 

Ji 

Once  we  have  estimated  J  the  optimum  choice  for  the  gab  is  /f  =  7/  where  7  is  a  real  scalar  and 
the  overbar  denotes  complex  conjugate.  This  control  will  use  the  least  amount  of  control  energy  for 
a  fixed  amount  of  stabilization. 


2.5.  Simulation  Results 

We  consider  a  nonlinear  plant  which  imdergoes  a  Hopf  bifurcation  to  a  limit  cycle  under  the  action 
of  a  Ibear  regulator.  We  study  the  effect  of  controller  gain  (for  a  fixed  phase)  on  the  amphtude  of 
the  hmit  cycle  oscillations.  We  also  perform  a  computational  experiment  for  fixed  gab  and  vary  the 
phase  of  the  controller. 

First  we  consider  the  “steady-state"  solutions  of  equation  (4.6).  The  term  steady-state  b  this 
context  refers  to  a  steady  Fourier  amplitude  or  constant  “energy”  £  =:iy  |*  of  the  synstem  Fbr  a 
pitchfork  bifurcation  this  is  a  steady-state  b  time.  For  a  Hopf  bifurcation  this  corresponds  to  a 


time-periodic  state  or  limit  cycle.  We  write  KG(iu>^°^)  =  Ke**  so  that  the  gain  of  the  compensator 
is  unity  at  the  critical  frequency  (|G(iu;(°))  [=  1).  When  the  system  is  in  a  lim  .t  cycle,  the  steady 
state  solution  of  equation  (4.16)  yields 
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\Y\= 


/ 


—  JrK  coe4>  +  JiKtxa(l> 

Sr 


(5.1) 


We  cannot  uniquely  determine  the  parameters  Jj  and  using  steady-state  data.  If  one 

normalizes  by  the  amplitMde  without  control,  jV  |nc=  the  normalized  amplitude, 

|yNy|/|y  Inc,  satisfies 


\y\= 


Jr 

AOi 


K  coe4  + 


^K^4> 


(5.2) 


That  is,  we  normalize  the  coefficients  by  the  growth  rate,  A^^,  which  cannot  be  detemuned  by 
steady-state  measurements.  The  normalized  amplitude  exhibits  the  foDowmg  dependence  on  the 
controller  gain  and  phase: 

|yi=  y/l  -  K{acoa(j>- P^4>)  (5.3) 


where  a  and  ^  depend  on  the  plant  parameters  via  equation  (4.6).  A  qualitative  prediction  from 
(5.3)  is  that  for  a  fixed  small  gain,  it  is  easier  to  reduce  limit  cycle  oscillations  than  it  is  to  enhance 
them  using  harmonic  forcing.  This  is  verified  by  simulations  in  this  section  and  by  fluid  experimenU 
in  the  next  section.  For  large  gains,  on  the  other  hand,  one  generally  destabilizes  another  mode  and 
this  relation  no  longer  holds. 

Example:  Van  der  Pol  OtcUlaior  with  a  Linear  Regulator  Simulations 

We  now  describe  the  nonlinear  plant  model,  the  Vwi  der  Pol  oscillator,  in  state  space  form. 

Xpi  —  -Xps  +  Bfiu 

Xp3  =  ~Xpi  -f  fiXpi  —  XpiXpj  -f  BpiU  (5,4) 

V  ~  Cp\Xpi  -f  CpiXpi 

The  linear  regulator  consist  of  a  gain,  a  band  pass  filter  and  a  lead  compensator  in  series  used 
as  a  linear  regulator,  that  is  G(s)  =  B{s)L{s).  The  band  pass  filter  has  the  transfer  function 


B{s) 


es 

s*  -f  cs  -f  1 


(5.5) 


where  the  peuameter  c  determines  the  bandwidth  of  the  filter.  The  lead 
function 


L{s)  = 


T»(s-n/r)’ 

(s  +  T)3 


compensator  has  the  transfer 


(5.6) 


The  maximum  lead  occurs  at  1  rad/ tec  and  its  value  b  LL[i)  =  2tan~*(i(7'-  l/T)).  By  changing 
the  sign  of  the  gain  we  can  obtain  phase  changes  from  -ISOdeg  to  -flSOdeg  by  varying  T.  Note 
that  |G|=  if  and  d>  ~  LG  =  LL. 

The  values  used  m  the  simulation  are  Bp\  =  Qj  =  0,  Bp^  =  =  1,  c  =  .5  and  p  =  .1. 

Changing  the  values  of  Bp  and  Cp  will  only  change  Jip  so  that  the  gain  and  phase  of  the  compensator 
would  have  to  be  adjusted  accordingly.  The  results  are  the  same,  scaled  in  gain  and  shifted  in  phase. 
For  these  particular  values  suppression  of  limit  cycle  oscillation  occurs  at  90  deg  phase  and  gain 
K  =  fi.  According  to  the  asymptotic  theory,  the  normalized  amplitude  should  be  given  by 


l>'l= 


1 - ain^ 


(5.7) 
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The  data  from  the  simulatioDB  are  depicted  as  circles  in  Figure  3  and  the  solid  curve  is  the  value 
predicted  by  the  Stuart-Landau  equation  for  K  =  =  If  one  increases  the  gain  K  above  fi  then 

the  limit,  cycle  amplitude  will  be  suppressed  for  a  range  of  phases  symmetric  about  90  deg.  As  me 
varies  the  gain  of  the  compensator  with  a  phase  of  ±90  deg  the  normalized  amplitude  behaves  as 
expected,  see  Figure  4. 

An  important  result  for  the  practitioner  to  be  aware  of  b  that  using  a  band  pass  filter  which  b 
too  narrow,  will  violate  the  assumptions  used  in  deriving  the  Stuart-Landau  equation.  Note  that 
G{i)  b  independent  of  the  value  of  c.  A  narrow  band  pass  filter  results  when  c  becomes  smaller. 
As  one  decreases  c  the  filter  poles  move  toward  the  imaginary  axb.  When  c  ~  then  the  filter 
poles  are  of  the  same  order  as  the  critical  poles  of  the  plant,  so  that  they  would  need  to  be  accounted 
for  m  the  multiple  scale  analysb,  leading  to  an  additional  amplitude  equation.  In  practice,  thb  b 
realbed  as  a  degradation  in  performance.  Thb  b  shown  for  the  Van  der  Pol  oscillator  m  Figure  5. 
As  the  band  pass  filter  becomes  more  narrow,  the  limit  cycle  amplitude  b  not  suppressed  as  well, 
and  b  no  longer  well  predicted  by  the  multiple  scales  theory  (5.7). 

We  now  consider  control  induced  transients  to  estimate  the  growth  rate  Ar-  Tb  simulate  the 
estimation  of  the  parameters  in  the  Stuart-Landau  equation  for  an  unknown  plant,  we  pretend  that 
we  do  not  know  the  form  of  the  governing  equation  (5.4).  We  pick  an  arbitrary  small  gain  and  four 
arbitrary  phases  for  the  compensator  G(i).  Here  we  use  K  —  fi  =  0.1  and  ^  =  60, 120,  240  and  300 
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120* 

O 

O 

O 

O 

Least  Squares 

Af^ 

.0063 
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Estimate 
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Multiple  Scales 

A^^ 

.0067 

.0067 

.0933 

.0933 

Computation 

Sr 

.1250 

.1250 

.1250 

•.1250 

Ibble  1 ;  Parameter  estimates  for  closed  loop  Stuart-Landau  equation 


degrees.  The  uncontrolled  system  exhibits  limit  cycle  oscillations  of  magnitude  0.6  for  p  =:  0.1.  At 
time  1  =  0  the  control  b  turned  on.  The  results  are  shown  in  Figure  6.  Using  equation  (4.15)  we 
determine  the  estimate  of  the  coefficients  for  the  closed  loop  Stuart-Landau  equation,  and  pr, 
given  in  Tbble  1. 

We  list  the  true  values,  and  gr,  computed  from  the  governing  equation  (5.4)  for  comparison. 
As  a  general  rule,  the  ratio  of  K^^/gr  b  quite  close  to  the  true  value  A^^/gr  since  the  estimates  for 
the  limit  cycle  amplitude  are  good.  The  growth  rates  are  more  difficult  to  obtain. 

Using  the  values  in  Ihble  1  we  can  use  equation  (4.18)  to  find  the  estimates  of  the  coefficients 
for  the  open  loop  Stuart-Landau  equation,  Ar  and  J,  as 

Ar  =  .0312 

J  = -.0021  -1.2921  (5.8) 

The  values  computed  from  the  governing  equation  are 

Ar  =  .05 

y  =  0.0  -  iO.5  (5.9) 

From  the  estimated  values  of  the  prirameters  (5.8)  the  compensator  that  should  give  marginal  sta- 
bihty  has 

/fG(i)=  j^  = -.0008  +  1.1067  (5.10) 

that  b  K  =  .1067  and  <p  =  90.409  deg  as  opposed  to  the  exact  values  of  A  =  .1  and  4>  =  90  deg 
Using  this  value  for  the  compensator  we  obtain  the  first  time  series  in  Figure  7  The  second  time 
series  b  for  the  same  phase  but  with  a  gain  K  =  .12. 
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Figure  7:  After  parameter  estimation  we  may  select  a  gain  and  phase  which  stabilises  the  plant. 


2.6.  Experimental  Results 

To  test  the  relevance  of  the  control  scheme  to  fluid  flow  oscillations  an  experimental  apparatus 
has  been  studied.  The  plant  we  desire  to  control  is  a  bested  two  dimensional  jet.  The  details  oi 
the  hot  jet  facility  and  some  of  its  uncontrolled  characteristics  are  discussed  in  [12].  It  has  been 
demonstrated  that  the  hot  jet  exhibits  self-excited  oscillations  when  the  ratio,  S,  of  the  jet  density 
to  the  ambient  fluid  density  is  below  a  uitical  value  5c  »  .95.  A  linear  stability  analysis  of  this 
flow  is  given  in  [12]  and  [2].  Above  this  ratio  the  flow  is  lammar  when  there  is  no  noise  present.  As 
the  density  ratio  is  decreased  the  steady  laminar  flow  undergoes  a  Eopf  bifurcation  to  a  limit  cycle. 
The  self-excited  nature  of  the  hot  jet  is  similar  to  that  of  the  von  Karman  vortex  street  in  the  wake 
behind  a  cylinder. 

ib  verify  that  the  hot  jet  exhibits  limit  cycle  oscillations  which  are  consistent  with  our  assump¬ 
tions,  we  measured  the  spectral  peak  for  several  different  temperatures.  The  data  presented  in 
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Figure  8  shows  experimental  data  acquired  on  three  different  days  (  o,  x,  and  +)  and  the  least 
squares  fit  of  the  data  to  the  equation  for  an  imperfect  bifurcation 

(^)|i'|-Piiy|"-P2  =  0  (6.11) 

To  implement  the  control  we  choose  acoustic  forcing  which  acts  like  a  point  actuator  at  the  lip 
of  the  jet  according  to  Bechert  [13, 14]  and  Crighton  [15]  who  have  shown  that  long  acoustic  waves 
can  be  converted  to  short  vorticity  waves  where  the  flow  has  a  strong  spatial  inhomogeneity.  The 
control  actuators  consist  of  two  100  Watt  woofers,  <«e  on  each  side  of  the  jet  approximately  1  meter 
downstream  of  the  jet  exit.  For  commercial  applications,  one  would  likely  use  some  form  of  vibrating 
plate  at  the  nozzle  exit  (for  an  example  of  such  an  actuator  see  [16]).  In  order  to  demonstrate  the 
feasibility  of  the  control  systems  considered  here,  the  acoustic  actuators  are  however  effective  and 
much  less  expensive. 

The  single,  point  sensor  is  a  hot  wire  probe.  The  anemometer  is  operated  in  constant  current 
mode  such  that  the  output  of  the  measurement  conditioning  circuit  is  proportional  to  the  fluctuating 
component  of  the  temperature,  J^(x,,t).  The  controller  consists  of  analog  filters  and  an  Intel  486 
based  personal  computer.  Rather  than  separately  estimate  the  dynamics  of  the  speakers,  the  time 
delay  associated  with  the  speed  of  sound  and  the  transfer  functions  of  various  filters,  all  the  dynamics 
are  lumped  together  and  considered  as  part  of  the  (unknown)  plant.  The  problem  is  then  reduced  to 
estimating  the  plant  (i.e.  the  coefficients  of  the  Stuart-Landau  equation)  and  specifying  the  control 
system  gain  and  phase  at  the  critical  frequency. 

The  sensor  was  located  at  xfH  —  1.4,  yJH  k.  0.5  and  2  =  0  which  is  upstream  of  the  initial  roll 
up  of  the  vortices.  A  probe  placed  in  the  shear  layer  at  this  location  has  a  distinct  spectral  peak  at 
about  100/f  2.  As  the  probe  is  moved  downstream  inside  the  shear  layer  one  enters  a  region  of  mode 
switching  in  which  the  Fourier  amplitudes  at  50if2  and  lOOHz  exchange  energy  in  time.  This  is 
believed  to  be  a  result  of  vortex  pairmg  which  is  not  completely  fixed  b  space.  Further  downstream 
the  vortex  has  paired  and  the  primary  feature  of  the  flow  is  a  large  peak  m  the  spectrum  at  SOHz, 
and  a  smaller  peak  at  100/f  2. 

We  now  consider  the  results  of  experiments  on  the  hot  jet  with  a  feedback  controller.  We  specify 
the  gab  and  phase  of  the  controller  at  the  critical  frequency  of  the  uncontrolled  jet.  The  data 
in  this  section  are  steady-state  b  the  sense  that  the  Fourier  amplitudes  of  the  critical  frequency 
(approximately  100/f  z)  are  steady,  that  is  the  measured  output  p(<)  is  time  periodic. 

Recall  that  the  mode)  for  the  flow’s  oscillations  is  the  Stuart-Landau  equation  with  a  control 
term  added 

^  =  K°^Y-g\Y\^Y-VJU  (6.12) 

where  U  =  Ke**Y  =  KG{iu^^^)Y  for  feedback  control.  Sbee  the  flow  b  spatially  dependent,  the 
coefficients  J  and  g  will  depend  on  the  sensor  location  x,.  Sbcc  results  from  the  real  part  of 
the  critical  eigenvalue,  it  must  be  bdependent  of  the  sensor  location. 

As  discussed  in  Section  2.5,  we  can  study  the  validity  of  the  Stuart-Landau  equation  as  a  model 
for  the  plant  by  examinbg  steady-state  data.  The  Fourier  amplitudes  ere  computed  from  the  hot 
wire  output  b  constant  current  mode  using  a  spectrum  analyzer.  A  distbet  spectral  p)eak  is  present 
for  a  oensity  ratio  of  S  <  Se  rs  .95.  At  a  density  ratio  of  S  —  .88  the  effect  of  controller  phase  has 
been  tested  and  compared  to  the  predictions  of  the  Stuart-Landau  model  b  Figure  9.  The  probe 
is  located  bside  the  shear  layer  and  the  amplitude  b  the  figure  is  normalised  by  the  open  loop 
amplitude.  The  solid  curve  is  a  least  squares  fit  of  the  data  to  the  Stuart-Lr  au  model 

l>'r=Pi +P3C06^-f  pssb^  (6.13) 

The  model  is  successful  at  predietbg  the  response  of  the  hot  jet  to  different  Ibear  regulators. 

The  effect  of  controller  gain  is  shown  b  Figure  10  for  the  same  plant  parameters  as  were  used 
to  study  the  effect  of  controller  phase  The  data  is  curve  fit  to  the  equation  (6.11)  for  an  imperfect 
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bifurcation  with  S  replaced  by  K  and  pi  and  P2  consistent  with  the  normalization  of  Y  to  3’- 
The  effect  of  random  disturbances  and  unmodeled  dynamics  prevent  the  oscillations  from  being 
completely  suppressed  for  any  value  of  the  controller  gain.  The  dashed  line  is  the  least-squares  fit 
to  the  imperfect  bifurcation  and  the  solid  line  is  the  perfect  bifurcation  (i.e.  P2  —  0). 

Another  prediction  of  the  multiple  scales  theory  for  small  gmn  controllers  is  that  the  eigenfunction 
for  the  least  stable  eigenvalue  is  the  same  with  control  as  for  the  uncontrolled  case.  That  is,  the 
spatial  structure  is  the  same  with  or  without  control.  This  is  verified  in  Figure  11.  The  figure  depicts 
a  Schlieren  inoage  with  and  without  control. 

In  order  to  estimate  the  growth  rate  Ar  in  the  Stuart-Landau  equation  we  must  perform  transient 
experiments  comparable  to  Section2.5.To  obtain  time  series  suitable  for  estimation  of  the  coefficients 
in  the  Stuart-Landau  equation  we  start  with  the  controller  turned  on  at  phase  which  causes  the 
maximum  amplification  in  Figure  9.  At  time  f  =  0  a  digital  switch  switches  the  phase  by  180  deg 
(i.e.  multiplies  the  control  signal  by  —1).  The  resulting  transient  is  shown  in  Figure  12.  The  first 
time  trace  is  raw  data  and  the  second  time  series  is  bandpass  filtered  and  displays  the  amplitude 
obtained  using  the  Hilbert  transform. 

Fifty  of  these  transients  were  obtained  and  averaged  to  obtain  an  average  for  the  amplitude 
|y(t)|  and  the  instantaneous  frequency  ^Ziy(t).  The  result  is  displayed  in  Figure  13.  Note  that  the 
frequency  shifts  and  becomes  somewhat  erratic  when  the  controller  switches  to  the  phase  for  reduced 
oscillation.  This  inhibits  our  ability  to  completely  suppress  the  flow  oscillations.  If  we  increase  the 
gain  further  we  find  that  a  different  mode  is  destabilized:  While  the  noode  at  lOOHz  is  suppressed, 
a  mode  at  150/fz  appears.  This  result  is  generic  and  has  also  been  demonstrated  for  the  wake  in 
[11- 

In  summary,  a  procedure  for  determining  the  effect  of  a  linear  regulator  on  bifurcating  solutions 
of  weakly  nonlinear  ordinary  differential  equations  has  been  developed.  Using  this  model  one  can 
specify  the  gain  and  phase  of  a  linear  controller  to  stabilize  the  system.  One  may  also  specify  a 
controller  to  modify  the  size  of  the  limit  cycle. 

Often  one  does  not  know  the  differential  equations  governing  the  system,  but  one  does  know  that 
the  system  exhibits  behavior  similar  to  a  Hopf  bifurcation.  The  Stuart-Landau  equation  can  then  be 
used  as  a  generic  model  for  the  system.  One  can  then  estimate  the  parameters  of  the  Stuart-Landau 
equation  using  experimental  data.  Different  controllers  can  be  used  to  provide  transient  data  for 
estimation  without  the  need  to  change  the  plant  parameters  for  their  generation.  (9,  17,  18]. 

Simulations  of  the  Van  der  Pol  oscillator  and  experiments  on  a  heated  two  dimensional  jet 
demonstrate  the  usefulness  of  the  Stuart-Landau  equation  for  control  and  estimation.  Work  is  still 
underway  to  develop  a  fully  adaptive  control  scheme  using  the  Stuart-Landau  equation  as  a  reference 
model. 


Figure  8;  Hopf  bifurcation  for  the  heated  two  dimensional  jet 
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3.  IMPLEMENTATION  OF  THE  CONTROL  IN  THE  WAKE  OF  A  CYLINDER  AT  LOW 
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3.2.  The  Doubly- Infinite  GL  Model  with  Control 

In  the  following  we  implement  the  controller,  developed  in  Section  2,  in 
the  doubly -infinite  Ginzburg-Landau  (GL)  model  of  Section  1.4.  which  we  are 
using  for  a  cylinder  wake  at  low  Reynolds  number.  Since  the  GL  model  retains 
the  streamwise  structure  of  global  modes,  the  locations  and  x  of  the 
single  sensor  and  of  the  single  actuator,  respectively,  have  to  oe  specified. 
Furthermore,  as  it  has  been  shown  in  Section  2,  the  controller  gain  and  phase 
are  only  relevant  at  the  frequency  of  the  critical  (most  amplified)  mode. 
Therefore  it  is  possible  to  use  a  simple  proportional  feedback  on  equation 
(4.9)  of  Section  1. 

Reverting  to  unsealed  physical  coordinates,  the  GL  equation  for  the 
amplitude  A  of  any  physical  quantity  is  rewritten  as 

a^A  -  L,  a^A  -  Lj  a^A  -  L^  a  +  N^  IAI^A  -  F(x,t)  .  (2.1) 

where,  using  the  notation  (4.10)  of  Section  1,  the  coefficients  are  defined  as 


Lo(^) 

—  -i[Wp(x)  +  5u)*’  + 

<))  ,  with 

a)o(x) 

-  +  (1/2)  - 

(x-x*")^  and 

ko(x) 

- 

•x^)  : 

L,(x) 

-  k/x)  ; 

^(x) 

-  (^1^/2)  : 

Ns(x) 

-  n^  +  n^  -  complex  constant 

with  n  >0 

r 

In  equation  (2.1)  we  have,  in  addition,  introduced  the  standard  nonlinear 
term,  which  can  bj  obtained  from  a  weakly  nonlinear  extension  of  the  analysis 
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of  Section  1  (work  in  progress  by  LeDizes  and  Monkewitz),  fnd  a  forcing  term 
F(x,t).  To  represent  the  controller  designed  in  Section  2,  the  forcing  term  is 
specified  as 

F(x,t)  -  g  exp(i7)  6(x-x^)  A(x^,t)  ,  (2.3) 

where  g  is  the  controller  gain  and  7  the  controller  phase.  We  note  that  the 
forcing  term  (2.3)  preserves  the  eigenvalue  character  of  the  problem.  At  this 
point  we  consider  specifically  the  wake  of  a  oblonp  limorph  (piezoceramic) 
cylinder.  By  applying  a  voltage  across  electrodes  .  top  and  bottom  of  the 
cylinder,  it  can  be  moved  in  the  direction  normal  to  the  oncjming  flow  and  is 
used  directly  as  an  actuator  in  our  setup.  Hence  x^-0  and  (2.3)  expresses  the 
fact  that  the  sensor  signal  is  amplified,  phase-shifted  and  then  fed  directly 
to  the  actuator  after  "mild"  bandpass  filtering.  The  latter  had  typical 
cutoffs  at  0.5  and  1.5  times  the  dominant  Karman  frequency,  which  was 
sufficient  to  implement  the  control  of  Section  2  because  the  dominant  spectral 
peak  in  the  low-Reynolds  number  wake  is  at  least  20dB  above  all  other  spectral 
features.  The  oblong  cylinder  of  thickness  D“0.69mm  ,  a  chord  of  T“1.68mm  and 
a  length  of  approximately  lOOD  is  shown  schematically  on  Figure  la,  together 
with  the  position  of  the  hot  wire  sensor  at  x/D-10  and  y/D-1.  The  cylinder  was 
mounted  across  the  nozzle  of  a  high  quality  jet  facility  as  shown  on  the 
photograph  of  Figure  lb. 


Figure  1.  a)  Schematic  of  the  experimental  arrangement  with:  C,  oblong 
cylinder;  HW,  hot  wire,  b)  View  of  the  oimorph  cylinder  mounted  on  the  nozzle 
of  the  jet  facility. 
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Next,  the  global  stability  characteristics  of  the  uncontrolled  oblong 
cylinder  wake  was  determined  by  transient  measurements  analogous  to  the  ones 
of  Raghu  &  Monkewitz  (1991)  in  a  hot  jet  (see  also  Schumm,  1991).  The  results 
are  shown  on  Figure  2  and  the  correlations  for  the  global  linear  growth  rate, 
the  linear  frequency  and  the  saturation  (limit  cycle)  frequency  are  given  in 
the  caption. 


Figure  2.  Experimentally  determined  global  stability  characteristics  of  the 
bimorph  cylinder  wake,  a)  Linear  growth  rate  a-^'Xr/v  with  curve  fit 
(7-0. 123(R-R^)  .  b)  O,  linear  frequency  with  fit 

f 27+0. 142 (R-Rq)  :  □,  saturation  frequency"  f  -f*  with  fit 

f  -9.27+0.177(R-r")  .  “ 
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From  the  results  of  Figure  2  and  an  estimate  of  the  length  of  the 
recirculation  region  in  the  near-wake  coupled  with  local  stability 
calculations,  the  parameters  (2.2)  have  been  estimated  as  follows: 


R  -  79.5 

C 

x*”  -  (1.183  -  0.031i) 

+  Sw^  -  (0.708  +  0.081i)  +  (0.0015+  O.Olli)x(R-R^) 

-  (1.452  -  0.844i)  +  (0.018  +  0.019i)x(R-R^) 
-(0.00175  +  0.04791) 

(0^  -  (0.107  -  0.065i) 

V"r  -  -1.8 


(2.4) 


The  ratio  is  a  universal  constant  for  a  given  system  and  is  directly 

related  to  the  nonlinear  frequency  shift.  The  value  of  n^,  on  the  other  hand, 
can  be  chosen  arbitrarily  and  is  related  to  the  normalization  of  the 
amplitude.  It  is  noted  that  the  above  parameters  have  not  yet  been  obtained 
from  the  solvability  conditions  developed  in  Section  1  due  to  a  lack  of 
precise  mean  flow  Information,  but  represent  an  a  priori  estimate  which  has 
not  been  adjusted  to  fit  the  experiments. 


With  (2.4)  the  linear  stability  of  time-harmonic  solutions  of  equation 
(2.1)  with  homogeneous  boundary  conditions  at  IxJ-*®  and  the  control  (2.3)  can 
be  investigated.  As  long  as  all  free  (uncontrolled)  global  modes  are  damped 
the  control  can  only  be  used  to  destabilize  the  system  for  some  phase  and  a 
gain  beyond  a  critical  gain.  The  situation  becomes  more  interesting,  when  the 
first  global  mode  (n-0  in  equ.  4.15  of  Section  1)  becomes  unstable.  In  this 
case  the  controller  with  nonzero  gain  is  needed  to  stabilize  the  fundamental 
Karman  mode  in  the  wake.  The  linear  stability  boundaries  for  the  wake 
parameters  (2.4)  and  supercritical  Reynolds  numbers  are  shown  on  Figure  3. 


Figure  3.  Stability  boundaries  predicted  from  the  GL  model  with  the  parameters 
given  by  (2.4) . 
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The  lower  stability  boundary  on  Figure  3  represents  the  MINIMUM  gain  g 
required  to  stabilize  the  Karman  mode  at  the  OPTIMUM  controller  phase.  The 
upper  boundary,  on  the  other  hand,  represents  the  MAXIMUM  gain  below  which  the 
system  can  be  stabilized  for  at  least  one  controller  phase.  In  other  words, 
beyond  the  upper  boundary  at  least  one  (controlled)  global  mode  is  unstable  at 
ANY  controller  phase.  A  typical  situation  on  the  upper  boundary  is  shown  on 
Figure  4:  Although  the  n-0  (Karman)  mode  can  be  stabilized  for  a  small 
interval  of  controller  phases  around  w,  the  range  of  phases  for  which  a 
"higher  mode"  is  destabilized  completely  overlaps  the  stabilization  interval 
of  the  Karman  mode.  What  becomes  clear  from  the  result  of  Figure  3  is  that  the 
range  of  Reynolds  numbers  for  which  our  controller  is  able  to  stabilize  the 
system  is  rather  limited.  The  reason  for  this  is  that  the  controller  is 
designed  to  suppress  only  a  single  mode,  namely  the  Karman  mode,  under  the 
assumption  that  all  other  modes  are  well  damped  (see  Section  2),  which  is 
clearly  no  longer  the  case  in  the  GL  system.  The  fact  that  the  controller 
still  "works"  close  to  the  bifurcation  is  explained  by  the  very  small  gain 
required  to  stabilize  the  weakly  amplified  Karman  mode . 


O  T  2r 


Figure  4.  Typical  situation  for  the  temporal  amplification  rate  o  versus  phase 

on  the  upper  stability  boundary  of  Figure  3.  - ,  amplification  rate  of 

uncontrolled  Karman  mode.  U,  regions  of  instability. 


The  situation  becomes  even  clearer  if  equation  (2.1)  with  the  nonlinear  cem 
and  the  controller  (2.3)  is  integrated  numerically  to  yield  the  limit  cycle 
amplitude  as  a  function  of  controller  gain  and  phase.  The  resulting  amplitudes 
for  a  supercritical  Reynolds  number  of  81  are  shown  on  Figure  5a  on  which  the 
triangular  region  of  stabilization  in  the  gain-phase  plane  is  clearly  visible 
Figure  5b  shows  the  corresponding  limit  cycle  frequency  which  illustrates  the 
mode  switching  on  the  stability  boundary. 
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Figure  5.  a)  Perspective  view  of  the  controlled  limit  cycle  amplitude  versus 
controller  gain  and  phase,  b)  Perspective  view  of  the  difference  between 
controlled  limit  cycle  frequency  and  uncontrolled  linear  frequency  versus 
controller  gain  and  phase.  R-81,  Xj-O,  x^-2  . 


3.3.  The  Experimental  Verification  of  the  Model  Predictions 


The  model  predictions  of  Section  3.2.  have  been  tested  in  the  cylinder 
wake  of  Figure  1.  First  it  was  confirmed  that  the  subcritical  system  could  be 
destabilized  by  the  controller.  This  is  illustrated  in  Figure  6  by  the 
transient  of  the  sensor  signal  after  closing  the  feedback  loop.  The  processing 
of  the  transient  shown  on  Figure  6  b&c  is  described  in  detail  in  Raghu  and 
Monkewitz  (1991)  (see  also  Schumm,  1991)  and  yields  the  global  linear  growth 
rate  as  well  as  the  linear  and  the  saturation  frequencies.  In  an  analogous 
manner  it  is  shown  on  Figure  7  that  the  unstable  Karman  vortex  street  can  be 
stabilized  by  our  controller.  In  this  latter  case  it  is  noted  that  the  problem 
of  extraneous  noise  has  become  much  more  severe  and  starts  to  push  the  data 
reduction  scheme  to  its  limits.  It  appears  that  this  problem  with  the  signal- 
to-noise  ratio  is  pervasive  in  all  self-excited  OPEN  flows. 
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Figure  6.  For  caption  see  next  page. 


Figure  6.  a)  Transient  from  subcritical  (stable)  state  to  a  limit  cycle  after 
closing  the  feedback  loop  (linear  scale  on  top  and  logarithmic  scale  on 
bottom).  R-79.2,  gain  g-0.01  .  b)  Instantaneous  growth  race  of  transient  6a 
versus  amplitude  squared  A^,  yielding  0-O.2  .  c)  Instantaneous  frequency  of 
transient  6a  versus  A^,  yielding  f,  -8.48  and  f  -8.62  . 
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Figure  a)  Transient  from  supercritical  (unstable)  state  to  a  noisy  stable 
state  after  closing  the  feedback  loop  (linear  scale  on  top  and  logarithmic 
scale  on  bottom).  R-83.8,  gain  g-0.003  .  b)  Instantaneous  growth  rate  of 
transient  7a  versus  amplitude  squared  A**,  yielding  <7“»1  .  c)  Instantaneous 
frequency  of  transient  7a  versus  A^,  yielding  fj^^j^-9.6  and  f^^^-10  . 


The  stability  properties  of  the  wake  have  been  systematically  explored 
at  several  supercritical  Reynolds  numbers.  The  procedure  was  first  to  find  the 
optimal  phase  for  suppression  of  the  Karman  mode,  and  then  to  increase  the 
gain  at  that  fixed  optimal  phase,  while  recording  the  spectrum  of  the  sensor 
signal  on  an  HP  spectrum  analyzer.  An  example  of  the  results  for  R-83.8  is 
shown  on  Figure  8.  The  bottom  part  of  the  figure  shows  the  magnitude  of  the 
dominant  spectral  peak  u*^  of  the  sensor  signal.  To  provide  a  direct 
physical  interpretation  of  gain,  it  is  defined  as 

g-v'rjj,j(cylinder)/u't^^(sensor)  ,  where  is  the  transverse  velocity 

amplitude  of  the  bimorph  cylinder  which  was  estimated  optically.  The 
bifurcation  values  g^  of  the  gain  were  determined  by  fitting  the  experimental 
data  with  the  steady-state  limit  cycle  amplitude  of  a  NOISY  system.  Denoting 
the  amplitude  of  the  sensor  signal  by  |A|^,  the  latter  is  given  by 

a|A|^  -  i/|A||  +  a  -  C  with  a  -  p(g-g^)  ,  (3.1) 

where  a  represents  the  external  noise.  Hence  we  determine  at  each  Reynolds 
number  the  parameters  /i,  v,  g^  and  q  by  least-square  fitting  to  the 
experimental  results.  The  location  of  the  first  bifurcation  at  gj^-1.3  10*’  has 
been  verified  by  transient  measurements  like  the  ones  shown  on  Figures  6  and 
7,  which  allowed  to  measure  the  linear  growth  rate  a  directly  as  shown  in  the 
top  part  of  Figure  8.  This  figure  clearly  demonstrates  the  existence  of  a  gain 
window  from  gj^  to  g2^,“8.1  10'^  where  another  global  mode  is  destabilized  by 
the  controller,  and  qualitatively  confirms  the  GL  model  predictions. 


g  X  10-' 


Figure  8.  Spectral  peak  of  streamwise  velocity  fluctuations  at  x/D-10,  y/D-1 
versus  gain  g-v'  ^(cylinder)/u'  (sensor).  O,  Karman  mode  at  320  Hz; 

- ,  equation  (3.5);  - ,  equ.  (Ta)  with  o-O;  0,  "higher  mode"  at  298  Hz. 

- ,  tqu.  (3.1). 

Furthermore,  the  triangular  shape  of  the  region  of  stability  in  the 
gain-phase  plane  (see  Figure  5)  has  been  verified  near  the  critical  Reynolds 
number.  The  result  is  plotted  on  Figure  9  where  the  stability  boundaries 
represent  the  destabilization  of  "higher"  modes  since  the  Karman  mode  is 
marginally  stable  for  the  experimental  conditions  shown. 
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Figure  9.  Experimental  stability  boundary  as  a  function  of  controller  gain  for 
different  sensor  locations  (phases).  \  is  the  Karman  wave  length.  R-79.3aR^. 


Finally,  all  the  eyptrimental  results  for  the  stability  boundaries  are 
compiled  in  the  last  icigure  10  which  has  to  be  compared  to  Figure  3  obtained 
from  the  GL  model.  Although  the  qualitative  agreement  is  excellent,  we  note 
that  the  experimental  Reynolds  number  range  from  79.5  to  about  90,  over  which 
wake  oscillations  can  be  suppressed,  is  approximately  twice  as  large  as  the 
theoretically  predicted  range  on  Figure  3.  This  is  due  to  the  fact  that  the 
parameters  (2. A)  have  been  estimated  in  a  rather  crude  manner  and  that  no 
attempt  has  been  made  to  fine-tune  them  for  a  quantitative  match. 


Figure  10.  Experimental  stability  boundary  of  the  controlled  wake  in  the  gain- 
Reynolds  number  plane.  O,  gain  g-v' (cylin(.-^r)/u' (sensor)  (left  gain  scale) 
at  which  Karman  mode  is  suppressed;  □,  galr  g  at  which  "higher  mode"  is 
destabiljzed;  ■,  amplifier  gain  (right  gair  scale)  at  which  "higher  mode"  is 
destabilized. 
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1.  Introduction 

The  stability  of  spatially  developing  free  shear  flows,  such  as  jets  end  wakes.  Is 
studied  under  the  assumption  that  the  basic  flow  is  two-diaensional  and  incoapressible, 
and  that  its  streamwise  Jevelopment,  characterized  for  Instance  by  its  width  ani,  is 
‘slow*  on  the  scale  of  a  typical  instability  wavelength  x,  that  is  c  ■  (a/a)x(<u/dx)  «  1. 
This  laeans  that  the  basic  flow  depends  only  on  the  transverse  coordinate  y  and  the 
rescaled  ‘slow*  streamwise  coordinate 

X  -  IX  .  (1) 

Furthermore,  the  non-parallelism  of  the  flow  is  assumed  to  be  auinly  due  to  pressure  or 
body  forces  such  that,  in  the  case  of  viscous  flows,  the  Reynolds  number  is  not  directly 
related  to  i.  Under  the  additional  essuaptions  of  an  Infinite  flow  domain  without 
internal  boundaries  and  locally  stnble  flow  far  up-  and  downstream,  the  line'.r  global 
modes  have  been  studied  by  Chomaz,  Huerre  b  Redekopp  [1990a},  Huerre  (  Honkewitz  [1990] 
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«nd  Monl;ewitz  [1990],  where  global  nodes  are  understood  to  be  tine-periodic  perturbations 
of  the  basic  flow  which  vanish  at  infinity  in  all  space  directions.  The  above  authors 
have  shown  formally  that  in  flows  without  solid  boundaries  the  global  nodes  are  ’driven* 
by  a  ’wave-maker’  centered  around  X,,  which  is  in  general  complex  and  defined  by 

du  dta 

—  (k..X.)  -  0  .  -  (k  X.)  -  0  :  «  -  «(k.)  .  (2) 

dk  dX 

In  words,  X,  is  a  saddle  point  of  the  local  absolute  frequency  u^(X)  (the  frequency  of  the 
node  with  zero  group  velocity)  obtained  from  linear  parallel  stability  analyses  of  the 
local  velocity  profile  at  each  X.  In  physical  terns  the  global  response  is  driven  by  the 
’local  oscillator’  at  X,  with  frequency  u,  which  has  the  ’most  compatible  neighbors’ 
oscillating  at  the  same  frequency  (to  linear  order  in  X-X,).  The  linear  analysis  (a  paper 
on  the  full  derivation  from  the  governing  equations  by  Honkewitz,  Huerre  t  Chomaz  is  in 
preparation)  shows  that,  far  from  X,,  global  node  shapes  can  be  described  by  WKB 
approxinati'ins  [see  for  instance  Crighton  t  Caster,  1976].  The  connection  of  the  IflCB 
solutions  across  X,,  which  is  a  second  order  turning  point  of  the  problem,  is  shown  to  be 
described  by  a  linearized  complex  Ginzburg-tandau  (LCGL)  equation,  valid  in  an  0(»‘'*) 
neighborhood  of  X,.  The  resulting  eigen-frequencies  of  the  low-order  global  nodes  are 
found  to  be  within  0(«)  of  which  means  that  the  flow  must  contain  an  interval  of 
absolute  instability  on  the  real  X-axis  in  order  to  become  self-excited,  i.e.  to  support 
tine-amplified  global  modes. 

In  this  paper,  we  extend  the  linear  global-node  analysis  to  the  weakly  nonlinear 
regime.  Since  marginal  instability  of  a  global  node,  which  is  a  prerequisite  for  a  weakly 
nonlinear  approach,  is  generally  reached  far  from  the  traditional  parallel-flow  stability 
boundary,  this  represents  .»  nontrivial  extension  of  Stuart’s  parallel  theory  [see  e.g. 
Stuart,  1971].  So  far,  nonlinear  effects  on  global  nodes  have  been  studied  by  adding  a 
cubic  nonlinearity  to  the  LCGL  equation  in  an  ad  hoc  fashion  [Chomaz,  Huerre  t  Redekopp, 
1990b].  More  recently,  in  an  effort  parallel  to  ours,  LeDizes  and  Huerre  [LeDizes,  1990] 
have  formally  analyzed  several  situations  with  different  relative  iii<)ortance  of 
nonparallel  and  nonlinear  effects.  In  the  following,  we  concentrate  on  the  case  in  which 
the  spatio-temporal  evolution  of  global  nodes  is  equally  affected  by  nonlinear  and 
nonparallel  effects  in  the  neighborhood  of  the  ’wave  maker’  at  X^. 

2.  TTie  Qtnzburg-Landau  Equation  Qovoming  the  Evolution  ot  Global  Moctee 

In  this  section  we  sketch  the  derivation  of  the  nonlinear  CGL  equation  in  the 
region  around  X,.  As  in  the  linear  analysis  [see  e.g.  H  »  M,  1990],  the  x-coordinate  and 
time  are  rescaled  according  to 
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f  -  ,)  .  (3.) 

T  -  <t  .  (3b) 

For  the  present  nonlineer  anelysis  the  reference  location  x,  1$  conveniently  replaced  by 
its  real  part  x,  ^  (see  equation  5).  The  slow  ti»e  T  characterizes  the  evolution  of  a 
localized  initial  impulse  into  a  linear  global  node.  If  we  wish  the  nonlinear  saturation 
of  a  global  mode  to  take  place  on  the  same  time  scale  (assuming  the  Hopf  bifurcation  to 
be  supercritical),  the  maximum  amplitude  of  the  disturbance  must  be  of  order  0(e‘'’)  and 
the  stream  function  may  be  expanded  in  powers  of 

•  -  ♦,(y.X)  +  +  o<«’)  .  (4) 

As  usual,  the  degree  of  supercriticality,  i.e.  the  temporal  growth  rate  j,  has  to 
be  of  order  0(().  The  new  feature  is  that  nonlinear  effects  ar»  concentrated  around  the 
location  of  maximum  amplitude  of  the  linear  global  mode.  Henci,,  one  has  to  distinguish 
between  a  maximum  amplitude  outside  an  0(c‘'’)  neighborhood  of  X,,  and  a  maximum  within. 
In  the  first  case  the  amplitude  at  X,  is  exponentially  small  and  the  ‘wave  maker*  region 
remains  linear.  Focussing  on  the  second  possibility,  we  show  that  it  leads  to  the  CGL 
equation  for  the  global-mode  amplitude.  To  ensure  that  the  maximum  amplitude  is  reached 
within  0(«*'*)  of  X,,  one  has  to  limit  the  spatial  growth  rate  ,  to  0(<*").  Finally,  for 
the  maximum  amplitude  to  be  of  the  same  order  as  |eJ(X,),  X,  must  be  within  0(<)  of  the 
real  X-axis.  In  summary,  we  have 


For  cs0{l),  the  mean  flow  is  expanded  around 


i*. 


*.(y.X)  -  ♦.(y.X, ,)  +  — (y.X, ,)  +  < - -(y.X,  )  ♦ 

ax  2  ax* 


(S) 


(6) 


and  the  amplitude  expansion  (4)  is  introduced  into  the  Euler  or  Navier-Stokes  equations. 
At  order  0(<‘'*)  the  linear,  parallel  stability  problem  for  the  mean  velocity  profile  at 
X,  ^  yields 

♦  -  A((.T)  ♦(y;X  )  E(x.t)  exp!-*,  ,(  ■*  n,  ,T)  -  c.e.  . 

(7) 

E(x,t)  -  expllk,  ,(x-x,  ,)  -  lu,  ,t)  . 

where  c.c.  stands  for  complex  conjugate.  As  there  is  only  a  slight  flow  divergence  within 
(sO(l),  the  transverse  structure  of  is  frozen  and  given  by  the  linear  eigenfunction 
♦(y..x, ,).  At  higher  order,  sol, ability  conditions  are  enforced,  noting  that  secular  terms 
have  to  be  suppressed  already  at  0(<)  because  of  the  nonparallel  basic  flow.  This  yields 
at  0(<*'*)  the  CGL  equation  with  j-th  degree  pclynomials  f^({)  as  coefficients. 

dA  a*A  dA  -2(*,  -  0,  ,T) 

—  +  —  +  l?,(0  A  +  1  •  |A|*A  -  F({.T)  (8) 

ai  a(*  af 
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The  forcing  function  F({,T)  has  been  added  here  for  later  reference.  The  linearized 
version  of  (8)  with  f-0  yields  a  spectrin  of  closely-spaced  global  eigen-frequencies 
-  w, ,  +  «(tn,  j  +  o’ +  tin”)  with  quantization  at  0(«)  1  The  associated  linear  global 
■odes  are  products  of  a  Gaussian  and  Henalte  polynoalals  of  degree  n  [H  ft  H,  1990]. 


Before  closing  this  discussion,  a  conent  on  the  solution  far  away  fron  X,  Is  In 
order.  For  |X-X,|>0(1)  the  global  aode  anplltude  Is  exponentially  sull  coapared  to  Its 
■axlwin.  There,  takes  the  fona  of  linear  UKB  approxliatlons  [C  ft  G,  1976] 


1  ■ 

-  «*'*/»(«)  A*(X,T)  ♦‘(yiX)  axp  -  k*(X)clX  - 


(9) 


where  e^'  and  are  the  subdoailnant  solutions  far  up-  and  downstream  of  X,,  respectively 
[H  ft  H,  1990].  The  amplitudes  a*(X,T)  are  governed  by  first-order  POE’s  and  have  to  be 
matched  to  the  spatio-temporal  amplitude  a({,T)  obtained  from  (8).  This  reveals  the  only 
‘blemish*  of  the  present  analysis  as  this  matching  cannot  be  carried  out  analytically.  It 
can  be  argued  though,  that  the  only  play  the  role  of  ‘passive  tails*. 


3.  Application  to  Proportional  Faadback  Control  of  Vortax  Shmdding  Bahind  a  Cylindar 

To  Illustrate  the  predictive  power  of  the  CGL  model  we  present  a  comparison  with 
feedback  experiments  In  the  cylinder  wake,  which  will  remain  qualitative  In  this  paper. 
As  In  most  practical  setups,  one  actuator  and  one  sensor,  typically  located  downstream  of 
the  actuator,  are  specified.  For  simplicity  we  assume  that  both  are  within  an  of 
X,,  that  the  sensor  or  probe  measures  directly  A(f-f^.T),  and  th't  the  actuator  acts  at  a 
point  and.  In  the  transverse  direction.  Just  forces  the  eigenfunction  ,).  Such  a 
proportional  or  constant-gain  feedback  Is  represented  In  equation  (8)  by  the  forcing  term 

F({.T)  -  g  mpIIt)  «({•(,)  .  (10) 

The  feedback  gain  g  and  the  phase  shift  y  Incorporate  In  practice  both  the  (coaq)1ex) 
electronic  gain  and  the  ‘receptivity*  of  the  flow  to  the  actuator.  It  Is  noted  In  passing 
that  this  model  applies  to  physical  setups  with  arbitrary  lead-  or  lag  compensators,  as 
only  their  gain  and  phase  shift  at  the  ‘ctrrier  frtquency'  ^  are  relevant.  This  model 
(equations  8  and  10)  In  Its  linear  form  'las  already  been  studied  by  Honkewitz  [1989]  and 
has  shown  Interesting  generic  behavior:  When  all  global  amides  are  linearly  stable  at  zero 
gain,  feedback  destabilizes  the  system  beyond  a  critical  gain,  which  Increases  as  the 
system  becomes  more  stable.  The  main  point  Is  that  It  Is  always  a  higher  mode  (n>0) 
and  not  the  least  stable  mode  that  Is  destabilized  by  feedback.  Therefore,  In  the 
situation  of  an  n>0  mode  which  Is  unstable  at  zero  gain,  one  finds  a  ‘gain  window*  where 
g  Is  sufficient  to  suppress  the  growth  of  the  n*0  mode  (at,  say,  the  optimum  phase  angle 
7),  but  Is  still  below  the  critical  value  at  which  higher  modes  are  destabilized.  This 
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window  narrows  as  the  degree  of  supercriticality  (e.g.  the  Reynolds  nimber)  is  increased 
and  control  is  typically  lost  not  far  beyond  critical  conditions. 

To  test  these  results,  experiaents  have  been  conducted  In  the  wake  of  a  cylinder,  a 
flow  which  satisfies  the  nodel  assumptions  quite  closely  [see  e.g.  Monkewitz,  1988].  A 
’Bimorph-transducer*  of  thickness  0-0.699n),  chord  T-1.68nn  and  a  length  of  approximately 
lOOD,  already  described  by  Berger  [1967]  [see  also  Berger  I  Schun,  1988],  served  as 
cylinder  (see  figure  1).  By  applying  a  voltage  between  electrodes  on  top  and  bottom,  the 
cylinder  was  oscillated  in  the  transverse  bending  mode  and  served  as  actuator.  The  probe 
was  a  constant  temperature  hotwire  located  at  x/0«10  and  y/0*l.  Its  signal  was 
linearized,  passed  through  an  tmplifier/phase  shifter  and  fed  to  the  Bimorph-transducer. 


I  f 

c) 


X 

■O' 


Figure  1.  Sch«Mtlc  of  tnt  oiportaonto)  trrontMnt  olth  C  ,  06)004  cy)indor  ,  Ml  ,  hot  «irt 


For  the  experiments  we  focussed  attention  on  the  suppression  of  Kanman  vortex 
shedding,  i.e.  on  the  supercritical  regime  with  Reynolds  number  Re>Re,,„  ,  where 
oblong  cylinder  [B  &  S,‘  1988].  In  the  course  of  the  study,  the  main 
qualitative  predictions  of  the  model  were  confirmed;  First,  we  found  at  slightly 
supercritical  Reynolds  numbers  a  gain  window  in  which  the  fundamental  Karman  mode  could 
be  suppressed.  The  result  for  Re*83.8  is  shown  on  figure  2.  The  bottom  part  of  the  figure 
shows  the  the  magnitude  of  the  main  peak  u',^  in  the  streamwise  velocity  spectrum,  as 
measured  by  the  hotwire  probe,  versus  gain.  To  provide  a  direct  physical  interpretation 
of  the  gain,  it  is  defined  as  g*v‘,^(cyl1nder)/u’^(probe),  where  v'^  is  the  transverse 
velocity  amplitude  of  the  Bimorph  transducer,  which  was  measured  optically.  The  phase 
shift  7  was  held  constant  at  its  optimal  value.  The  bifurcation  values  g^,,,  were 
determined  by  fitting  the  data  with  the  steady-state  limit-cycle  amplitude  |A|  at 
obtained  from  (8)  and  defined  by 

6|A|  •  *<|A|’  *  a  -  0  ;  6  -  •  (ID 

To  account  for  the  external  noise  (freestream  turbulence)  that  keeps  exciting  the  Karman 
mode  beyond  its  neutral  point,  a  forcing  term  o  has  been  added  to  the  other  fitting 
parameters  p,  v  and  in  (11).  The  location  of  the  first  bifurcation  at  B,»1.3  10'* 
has  in  addition  been  verified  by  transient  experiments  in  which  the  global  linear  growth 
or  decay  rate  a  was  determined  directly.  The  figure  clearly  shows  the  gain  window 
extending  from  g,  to  g^'B.l  10'’  where  another  global  mode  is  destabilized,  as 
qualitatively  predicted  by  the  atodel.  In  addition,  we  also  confirmed  the  second  mam 
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Figure  2.  sptctrti  p««k  of  itraMiic  vtiocitr  <t 

y/e«l  vtrius  gitn  g«v*(C]r1tnder)/u'(prote)  O  ,  KirMrt  BO<te 

•t  3J0  Hj.  —  .  Pgu  (11), - ,  tgu  (11)  with  o»0 

«  ,  “ftlghir  •od*"  *t  J9J  Ki  j  01) 


prediction  that  control  Is  lost 
relatively  close  to  ,  by  showing 
that  suppression  of  vortex  shedding  In 
our  setup  Is  only  possible  between 

4.  CondutlotM  and  Acknowiadgamtnts 

This  study  deMonstrates  that  a 
relatively  slaple  weakly  nonlinear 
iK>del  of  wake  oscillations  Is  capable 
of  capturing  the  rather  complex 
dynamics  resulting  from  the  addition  of 
feedback.  The  next  step  will  be  a 
quantitative  comparison  between  the 
model,  with  coefficients  either 
determined  from  stability  calculations 
or  iMasured,  and  the  present 
experiments.  Further  ahead,  the  success 
of  Ffowcs  Williams  and  Zhao  [1989]  In 
suppressing  Kannan  vortex  shedding  at 
eight  times  Rej,„.  where  the  situation 
Is  greatly  complicated  by  turbulence, 
awaits  explanation. 
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Control  of  Hopf  Bifurcations 
Edward  A.  Henrich,  D.  L.  Mingori  and  P.  A.  Monkewitz 
University  of  California,  Los  Angeles 


1  Introduction 

Many  systems  in  nature  exhibit  self-excited  oscillations  or  limit  cycle  behavior.  The 
bifurcation  from  a  steady-state  to  a  limit  cycle  is  generally  the  first  step  in  the  tran¬ 
sition  to  chaos  or  turbulence.  For  a  given  physical  system  there  are  often  design 
constraints  which  specify  a  hxed  region  of  parameter  space  in  which  one  would  like 
to  operate.  If  this  particular  choice  of  parameters  exhibits  some  undesirable  dynamic 
behavior  (e.g.  limit  cycles  or  chaos)  then  the  goal  of  the  control  engineer  is  to  sup¬ 
press  this  behavior  while  still  operating  in  the  same  parameter  range.  A  distinction  is 
made  between  a  ‘control  parameter’  or  ‘bifurcation  parameter’  which  is  the  parameter 
one  uses  to  study  successive  bifurcations  and  a  ‘control  force’  which  is  a  prescribed 
function  of  time  used  to  modify  the  system’s  output.  As  an  exzimple,  consider  that 
one  might  desire  to  suppress  the  Von  Karman  vortex  street  in  the  wake  behind  a 
cylinder.  The  control  parameter  is  the  Reynolds  number  which  may  be  fixed  at  some 
supercritical  value,  where  at  the  desired  operating  speed  the  wake  is  no  longer  laim- 
inar.  The  control  force  could  be  an  induced  sound  field  designed  to  eliminate  flow 
oscillations.  This  method  uses  the  control  force  to  stabilize  the  system  in  contrast  to 
methods  which  modify  the  control  parameter  until  the  state  behaves  as  desired  (1). 

2  Control  System  Analysis 

For  a  physical  system  with  a  finite  number  of  degrees  of  freedom,  the  dynamics  may 
be  described  by  a  set  of  first  order  ordinary  differential  equations.  Written  in  vector 
form 


m) 

di 


-i  =  F(X(0;p)  +  Bu{0 

(2.1) 

y  =  CX-f  D 

(2.2) 

X  €  il"  and  p,ti,y  t  R 

(2.3) 

where  X  is  the  system  state,  F  is  the  vector  field,  u  is  the  control  force,  y  is  the 
measured  output  and  /i  is  a  control  parameter  in  the  problem.  Equation  (2.1)  is 
presumed  to  have  an  equilibrium  solution  Xo(p)  which  satis.'.es  F(Xo(p);p)  =  0.  For 
p  >  0  we  assume  the  system  undergoes  a  Hopf  bifurcation  to  a  limit  cycle.  The 
goal  of  the  control  system  is  to  return  the  system  to  its  equilibrium  state,  which  is 
unstable  for  p  >  0. 


1 


Lbear  CompcDMtor 
i.(,)=lplG(.)e(t) 


•(0 


Nonliiiear  Plant 
X  =F(X;»i)  +  B* 
»  =  CX  +B 


m 


Figure  1:  Block  diagram  of  the  control  system 

One  method  for  controlling  a  set  of  ordinary  differential  equations  which  depend 
on  a  control  parameter  p  is  to  let  the  value  of  the  control  parameter  be  governed  by 
a  differential  equation.  The  time  derivative  ft  depends  on  the  state  of  the  system  and 
a  gain  specified  by  the  designer.  If  the  gain  can  be  chosen  such  that  p  drops  below 
the  critical  value,  then  the  n  +  1  dimensional  system  will  be  stable.  This  procedure 
is  described  in  [1]. 

The  goal  of  this  research  is  io  control  the  system  by  external  forcing  when  the 
control  parameter  is  fixed  by  some  other  design  constraint.  The  control  system  which 
is  studied  is  the  linear  regulator  as  depicted  in  Figure  1.  The  differential  equations 
governing  the  control  system  are  such  that  the  Laplace  transform  of  the  impulse 
response  of  the  linear  compensator  is  fiG{s). 

The  action  of  the  control  system  is  studied  by  using  a  multiple  scale  perturbation 
technique  [2,  3).  We  consider  the  case  where  the  real  part  of  the  critical  eigenvalue 
is  O(e^).  In  this  case  it  is  natural  to  introduce  a  second  time  scale  t  =  c’f,  where 
p  =  c*sgn(p).  We  treat  the  two  time  scales  as  independent  variables  so  by  the  chain 
rule  the  derivative  becomes 

d  d  .d 

We  use  this  formula  and  expand  (2.1)  in  a  Taylor  series  in  X  -  Xo  and  p.  Upon 
equating  terms  multiplied  by  like  powers  of  <  we  obtain  a  sequence  of  inhomogeneous 
partial  differential  equations.  The  elimination  of  secular  terms  leads  to  a  solvability 
condition  which  governs  the  slow  time  evolution  of  the  amplitude  of  the  limit  cycle 
oscillations.  The  solution  on  the  center  manifold  is 


X(0  ~  Xo  +  c/l(T)c^‘‘r  +  £A(T)e-^‘‘r 


(2.5) 


where  is  the  imaginary  part  of  the  critical  eigenvalue  and  r  is  the  corresponding 
right  eigenvector.  The  solvability  condition  can  be  written  as  an  amplitude  equation 


AA 

—  =  (sgn(p)A  -  \BCt  (?(tu>c))24  '-g\A^  A 

aT 


(2.6) 
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where  A  is  determined  by  the  real  part  of  the  eigenvalue  of  the  Jacobian  of  F,  g  is  de¬ 
termined  by  the  nonlinearities,  and  the  receptivity  IBCr  is  determined  by  the  manner 
in  which  the  control  force  couples  to  the  measured  output  (1  is  the  left  eigenvector 
of  the  marginal  mode).  In  fluid  mechanics  this  type  of  equation  is  known  as  iLc 
Stuart-Landau  equation. 

An  important  thing  to  note  in  equation  (2.6)  is  that  the  stability  of  the  fixed 
point  Xo  is  governed  by  the  control  system  only  near  the  critical  frequency  Ug.  In 
pau’ticular,  the  entire  class  of  control  systems  with  the  same  gain  and  phase  at  the 
critical  frequency  are  equivalent.  The  particulzir  gain  and  phase  required  of  the  control 
system  at  the  critical  frequency  is  determined  by  the  stability  of  the  uncontrolled  plamt 
(i.e.  A)  and  its  receptivity  from  measured  output  to  control  force  (i.e.  lJ5Cr). 

3  Estimation  and  Control 

The  results  of  the  previous  section  demonstrate  the  effect  of  a  linear  regulator  on 
a  nonlinear  plant  which  undergoes  a  Hopf  bifurcation  to  a  limit  cycle.  In  deriving 
the  result,  we  assumed  that  the  system  equations  (2.1)  were  known.  If  the  system 
equations  were  known,  we  would  likely  use  one  of  the  standard  control  design  methods 
for  nonlinear  plants  (e.g.,  describing  functions  or  Lyapunov’s  direct  method).  The 
more  interesting  situation  is  when  the  model  for  the  system  dynamics  is  unknown, 
but  what  is  known  is  that  the  system  exhibits  limit  cycle  behavior.  In  this  case, 
it  is  useful  to  think  of  the  Stuart-Landau  equation  (2.6)  as  a  generic  model  for  the 
dynamics  which  are  asymptotic  to  the  center  manifold.  In  many  physical  applications 
the  measurement  consists  of  a  time  history  of  a  single  variable..  In  these  cases,  the 
state  X{i),  A{r)  and  the  eigenvector  r  may  all  be  unknown.  This  difficulty  is  easily 
overcome.  Defining  Y{{)  =  2£A(T)C’r,  the  leading  order  asymptotic  expression  for 
the  output  is 

y(0  =  iZe[>'(Oe’"'‘j  +  CXo  +  Z)  (3.1) 

(see  (2.2)  and  (2.5)).  Returning  to  the  physical  variables,  one  obtains  an  amplitude 
equation  which  is  useful  for  describing  the  measured  output  of  an  experiment  in  terms 
of  a  dynamical  system  on  the  center  manifold 

^  =  /<(sgn(p)A  -  \BCv  -  s  15'  I'  r  (3.2) 

or 

^^  =  (A-JG(iu..))y-5|y|M'  (3.3) 

where  A,  J  and  g  are  unknown  complex  constcints. 

To  estimate  these  constants  from  experimental  data  one  needs  a  method  for  ob¬ 
taining  y(t).  One  way  to  obtain  the  amplitude  signal  y(t)  would  be  to  create  a  90 
degree  phase  shifted  version  of  the  output  y(/)  as  y.  The  modulus  of  the  complex 
signal  y(i)  =  y(t)  -f  !y(<)  is  |y(0 1=1  ^^(0 1-  define  an  analytic  signal  as  a  complex 
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signal  of  a  real  variable  whose  real  and  imaginary  parts  are  Hilbert  pairs  [4].  This 
analytic  signal  is  just  an  extension  of  the  ‘rotating  vector’  or  ‘phasor’  used  in  circuit 
analysis  and  systems  engineering.  The  Hilbert  transform  is 

i  =  H{y]  =  ip.V.  f°°  =  1 .  v(<)  (3.4) 

TT  J—oo  I  ~  T  Til 

where  *  denotes  convolution  in  time.  Defined  this  way  the  Hilbert  transform  of  a  real 
time  signal  is  a  real  time  signal.  Taking  the  Fourier  transform  of  (3.4)  we  have 

}  =  ~i  8gn(u;).F{y}  (3.5) 

which  shows  that  the  Hilbert  transform  is  indeed  a  90  degree  phase  shifter.  From 
(3.5)  we  find  that  the  Fourier  transform  of  the  analytic  signal  is 

f  2T{y]  u)  >  0 

^■(j}  =  Hv)  =  0  (3.6) 

[  0  u;  <  0 

Hence  we  can  obtain  the  analytic  signal  y{i)  by  using  (3.6)  and  an  inverse  Fourier 
transform 

y{t)  =  +  sgn(a;)).?‘{y} }  (3.7) 

We  use  the  analytic  signal  to  find  the  ‘envelope’  of  the  output  and  its  ‘instantaneous 
phase  or  frequency’.  The  envelope  is  just  |y  |  and  the  instantaneous  phase  is  Ly.  So 
we  can  identify  our  complex  amplitude  as 

l'(l)=|i/|£'>  (3.8) 

where  4>  ly  -  The  Hilbert  transform  represents  a  means  of  obtaining  the 
demodulated  output  Y{t)  from  the  output  y{i)  off-line.  Equation  (3.3)  is  solvable 
for  Y{t)  as  an  analytic  function  in  terms  of  the  coefficients  A,  J  and  g  so  there  are 
several  means  to  estimate  the  parameters  of  the  Stuart-Landau  equation  once  we  have 
obtained  y(f)  from  y{t)  using  the  Hilbert  transform. 

The  estimation  and  control  scheme  may  be  demonstrated  using  the  Van  der  Pol 
oscillator  as  a  model  plant.  The  plant  is 

y  +  (y’-/^)y  +  y  ==  (3-9) 

The  plant  equation  cind  the  differential  equations  for  the  linear  regulator  are  inte¬ 
grated  using  a  fixed  step  fourth  order  Runge-Kutta  algorithm.  We  use  several  test 
values  of  G{iuc)  to  identify  the  system  parameters  A,  J  and  g.  Once  the  parameters 
have  been  estimated,  we  use  this  model  of  the  dynamics  on  the  center  manifold  to 
select  the  optimum  values  for  the  controller  gain  and  phase.  The  results  of  the  final 
controller  are  depicted  in  Figure  2.  The  solid  curve  is  the  measured  output  y{t)  and 
the  dashed  curve  is  the  control  force  u{t). 
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Figure  2:  Stabilization  of  limit  cycle  oscillations 

4  Conclusions 

A  procedure  for  finding  a  model  for  data  which  lies  on  the  center  manifold  has  been 
developed.  Using  this  model  one  can  specify  the  gain  and  phase  of  a  linear  controller 
to  stabilize  the  system.  One  may  also  specify  a  controller  to  modify  the  size  of  the 
limit  cycle. 

Previous  work  with  wakes  and  jets  have  shown  the  Stuart-Landau  equation  to 
be  a  good  model  for  the  dynamics  on  the  center  manifold  [5].  We  are  currently 
investigating  the  results  of  this  paper  using  a  heated  two  dimensional  jet  experiment. 

The  research  described  in  this  article  was  performed  under  AFOSR  contract  89- 
0421. 
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